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ABSTRACT 

These  pages  present  a  study  of  gaseous  diffusion  into  and 
turbulent  entrainment  from  steady  as  well  as  nonsteady  fully 
cavitating  flows  and  characterize  the  nonsteady  cavity  as  a  dynamic 
element  in  a  water  tunnel  test  section.  An  approximate  theoretical 
analysis  is  done  for  gaseous  diffusion  across  steady  axisymmetric  ' 
cavities.  A  comparison  with  experimental  results  reveals  that  the 
present  approximate  analysis  is  suited  for  better  estimates  of 
gaseous  diffusion  across  the  gas-liquid  interface  of  axisymmetric 
cavities  than  past  analyses. 

Water-tunnel  experiments  for  two-dimensional  flows  were 
conducted  on  cavities  behind  a  stationary  and  oscillating  "flat 
plate"  (wedge)  hydrofoil.  It  is  found  that  the  steady  cavity 
pressures  were  significantly  higher  than  the  vapor  pressure  of  water 
and  larger  cavities  tend  to  maintain  higher  cavity  pressures. 
Unsteady  cavity  pressures  were  found  to  be  significant  only  at  low 
frequencies  of  profile  oscillation.  Downstream  and  upstream  of  the 
cavity,  the  oscillating  cavity  was  observed  to  behave  like  a  dipole 
source  at  low  frequencies  of  oscillation  and  like  a  quadrupole 
source  at  high  frequencies  of  foil  oscillation.  An  analysis  of  the 
oscillating  cavity  lengths  revealed  that  for  sinusoidal  motions  of 
the  foil,  the  first  two  harmonics  of  the  cavity  motions  are 
significant. 

A  mathematical  model  is  then  developed  to  estimate  the 
nonsteady  gaseous  entrainment  behind  two-dimensional  steady  and 
oscillating  cavities.  This  model  takes  into  account  the  varying  gas 
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pressure  as  well  as  varying  cavity  lengths*  It  is  observed  that 
this  model  has  the  proper  steady  state  limit.  Steady  and  nonsteady 
gas  entrainment  rates  were  estimated  for  the  two-dimensional 
cavity.  It  is  found  that  gas  pressure  within  the  cavity  affects 
the  steady  entrainment  behind  the  cavity.  The  phase  angle  between 
the  foil  motion  and  cavity  motion  is  determined  by  a  physically 
driven  argument  that  is  seen  to  be  consistent  with  experimental 
observations.  The  nonsteady  gaseous  entrainment  behind  the  cavity 
is  seen  to  depend  only  on  the  cavity  volume  fluctuations  and  cavity 
gas  pressure  fluctuations. 
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CHAPTER  1 
INTRODUCTION 

1.1  Preliminary 

In  many  liquid  flow  systems,  some  form  of  unsteady  cavity  flow 
exists  at  one  or  more  locations  in  the  circuit.  Such  unsteady  cavity 
flow  may  arise  intentionally  as  in  the  case  of  water  tunnel 
experiments  on  nonsteady  cavity  flows  or  unintentionally  as  in  the 
case  of  high  speed,  high  performance  turbomachines  in  various  liquid 
flow  systems  such  as  nuclear  reactors  and  liquid  rocket  engines*  In 
any  case,  the  presence  of  an  unsteady  cavity  within  a  liquid  might 
create  disturbances  that  could  tend  to  make  the  system  unstable.  The 
propagation  of  such  disturbances  becomes  even  more  significant  when 
small  gas  bubbles  are  present  in  the  liquid,  wherein  the  acoustic 
celerity  of  the  gas-liquid  mixture  is  greatly  reduced.  A  classic 
example  is  the  "POGO"  instability  observed  in  liquid  rocket  engines 
(1,2).  In  closed  circuit  water  tunnel  operation  involving  unsteady 
cavity  flow  experimentation,  mass  oscillations  of  the  entire  circuit 
and  "breathing  of  the  tunnel"  have  been  reported  (3-5).  In  order  to 
obtain  valid  nonsteady  force  measurements  involving  cavity  flow  in 
closed  circuit  water  tunnels,  the  nature  of  the  unsteady  cavity  as  a 
disturbance  source  and  the  dynamics  of  the  water  tunnel  itself  must  be 
clearly  understood. 

It  is  the  purpose  of  this  thesis  to  examine  the  physical  aspects 
(such  as  gaseous  diffusion  into  the  cavity  and  its  turbulent 
entrainment  into  the  cavity  wake)  that  govern  the  nature  of  unsteady 
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cavities  so  that  system  responses  to  unsteady  cavity  flows  may  be  more 
accurately  predicted  and  related  to  force  measurements  on  oscillating 
struts  and  foils.  This  latter  aspect,  involving  the  dynamics  of  any 
particular  facility,  is  outside  the  scope  of  this  study,  however. 


1.2  Cavity  Aspects 

Whenever  the  local  pressure  in  a  region  within  the  liquid  falls 
below  the  vapor  pressure  of  the  liquid  at  the  ambient  temperature,  the 
liquid  evaporates  and  a  cavity  is  formed.  The  size  and  shape  of  the 
cavity  is  governed  by  the  cavity  pressure  and  the  shape  and  attitude 
of  the  body  on  which  the  cavity  is  formed.  The  cavity  is  usually 
characterized  by  a  quantity  known  as  the  "cavitation  number"  which  is 
defined  as: 


where 


a 


2  • 
l/2p  tr 

wm  m 


(l.D 


a  *  cavitation  number, 

pa  *  free  stream  static  pressure, 

pc  -  cavity  pressure, 

Pa  ■  free  stream  density  of  the  liquid, 
Ua  *  free  stream  velocity. 


r  In  a  purely  vaporous  cavity,  the  cavity  pressure  is  equal  to  the 

vapor  pressure  of  the  liquid  at  the  ambient  temperature  of  the  liquid. 

I  Gadd  and  Grant  (6),  among  others,  have  demonstrated  that  the  measured 

i  cavity  pressure,  p  ,  within  fully  developed  natural  cavities  exceeds 

c 

j  the  liquid  vapor  pressure,  py,  because  of  the  unavoidable  presence  of 


gases  dissolved  in  the  liquid.  These  gases  diffuse  through  the  cavity 
interface  and  set  up  a  partial  pressure  of  gas  p  ,  in  the  cavity  so 
that 

pc  “  Pv  +  PG  *  (1*2) 

Thus,  the  presence  of  dissolved  gases  in  the  liauid  affect  the 
cavitation  number  a  and  the  behavior  of  the  cavity  itself  to  a 
great  extent. 

The  influence  of  dissolved  gases  on  the  behavior  of  single 
bubble  cavities,  viz.,  growth  and  collapse,  has  been  the  subject 
of  many  previous  investigations.  Epstein  and  Plesset  (7)  studied 
the  stability  of  gas  bubbles  in  liquid  gas  solutions.  Parkin  and 
Kermeen  (8)  studied  steady  gaseous  molecular  diffusion  into  a 
microbubble.  They  Investigated  the  influence  of  gaseous  diffusion 
on  the  growth  of  small  bubbles  in  a  flow  field  containing 
dissolved  gases.  In  cavity  flows,  when  the  boundary  layer  on  the 
cavity  surface  is  turbulent,  Brennen  (9)  has  shown  that  turbulent 
gaseous  diffusion  into  the  cavity  may  be  more  significant  than 
molecular  diffusion.  Brennen  (9)  also  extended  his  two-dimensional 
theory  to  include  the  gaseous  diffusion  across  an  axlsymmetric 
cavity  by  wrapping  the  two-dimensional  solution  around  the  cavity. 

Several  experimental  investigations  have  been  devoted  to  the 
understanding  of  gaseous  diffusion  across  steady  cavities  (10-13). 

A  majority  of  these  studies  are  devoted  to  the  understanding  of 
axisymmetric  cavities  because  of  frequent  occurrence  of  such  shaoes 
in  marine  applications.  Swanson  and  O'Neill  (10)  investigated  the 
stability  aspects  of  ventilated  cavities  behind  circular  discs  in 
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|  steady  flow.  Cox  and  Clayden  (11)  studied  the  air  entrainment  behind 

I 

cavities.  Billet  and  Weir  (13)  measured  the  entrainment  rates 
behind  ventilated  axisymmetric  cavities. 

1 

The  unsteadiness  of  a  cavity  arise  in  two  distinct  forms.  The 
first  one  is  the  auto  oscillation  of  steady  cavities  (particularly  the 
ventilated  cavities)  that  is  prevalent  at  certain  cavitation  numbers. 
This  phenomenon  has  been  studied  in  great  detail  by  Song  (14)  and 
[  Sllbeman  and  Song  (15).  The  second  form  is  the  forced  oscillation  of 

^  the  cavity  due  to  the  oscillation  of  the  body  on  which  the  cavity  is 

formed,  l^pical  examples  are  a  cavity  behind  an  oscillating  hydrofoil 
and  a  cavity  on  a  propeller  blade.  The  oscillations  of  the  cavity  may 
be  classified  into  two  distinct  forms  of  motion  (16):  First,  the 
motion  of  the  cavity  wherein  the  volume  of  the  cavity  changes 
significantly;  second,  the  motion  of  the  cavity  wherein  the  shape  of 
the  cavity  changes  significantly  but  the  cavity  volume  remains 
essentially  constant.  As  the  cavity  volume  changes,  the  partial 
pressure  of  gas  inside  the  cavity  changes.  Since  the  time  scale  of  the 
process  of  vapor  release  and  absorption  is  very  much  smaller  than  the 
period  of  cavity  oscillations  typically  encountered,  the  vapor  pressure 
inside  the  cavity  essentially  remains  constant.  Thus  the  cavity 
pressure  within  an  oscillating  cavity  may  be  written  as 

Pc(e>  ■  Pv  +  PG<0  »  d.3) 

where  p^(t)  is  the  cavity  gas  pressure.  Then  the  cavitation  number 
o(t)  of  an  oscillating  cavity  is  given  by 
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PG(t) 

o(t)  -  on - t  ,  (1.4) 

1/2P.U. 

where  Is  the  cavitation  number  based  on  the  vapor  pressure  of  the 
liquid.  Jiang  and  Leehey  (17)  have  measured  the  cavity  pressure 
within  an  oscillating  cavity  behind  an  elliptical  wing  oscillating 
in  pitch.  They  found  that  the  unsteady  cavity  pressure  depends  on 
the  frequency  of  oscillation  and  that  the  unsteady  pressures  are 
significant  at  low  frequencies  of  oscillation. 

One  of  the  earliest  models  used  in  characterizing  the  unsteady 
cavity  as  a  source  of  disturbances  in  fluid  flow  system  was  in 
connection  with  the  analysis  of  "POGO"  instability  in  liquid  rocket 
engines.  By  defining  a  quantity  called  "cavitation  compliance"  which 
relates  the  mass  flow  fluctuations  and  the  corresponding  pressure 
fluctuations,  Rubin  (2)  was  able  to  model  the  oscillating  cavitatlng 
bubble  as  a  simple  spring  mass  system.  Brennen  and  Acosta  (18)  showed 
that  Rubin's  model  for  a  cavitatlng  pump  inducer  is  rather  simple. 

The  results  of  their  experimental  investigation  revealed  appreciable 
departure  from  Rubin's  quasi-steady  model  due  to  dynamic  effects  at 
high  frequency.  Huse  (16),  in  an  attempt  to  predict  cavitation  damage 
on  the  hull  of  a  ship  due  to  a  cavitatlng  propeller,  represented  the 
unsteady  cavity  at  the  tip  of  the  propeller  as  a  combination  of 

* 

) 

acoustic  monopole  and  dipole  sources.  His  analytical  results  compared 
well  with  the  experiments. 


1.3  Scope  of  Present  Investigation 


The  gas  pressure  within  the  cavity  is  a  quantity  of  fundamental 
importance  in  steady  and  nonsteady  cavity  flows.  The  gas  pressure 
within  the  cavity  dictates  the  amount  of  gas  diffusion  across  the 
cavity  surface.  Hence  we  seek  to  describe  the  diffusion  phenomeon 
analytically  and  to  compare  the  results  with  the  best  available 
experimental  data.  To  our  knowledge,  the  data  of  Billet  and  Weir 
(13)  are  the  best  available  and  they  pertain  to  axlsyrametric  flows. 

In  Chapter  2  we  first  demonstrate  using  a  heuristic  model,  the 
lack  of  correlation  between  existing  mathematical  models  for  gaseous 
diffusion  across  axisymmetric  cavities  and  corresponding  experimental 
results.  Then,  we  use  the  convective  diffusion  theory  proposed  by 
Parkin  (8)  with  the  molecular  diffusivity  being  replaced  by  a  turbulent 
diffusivity  to  predict  the  gaseous  diffusion  across  an  axisymmetric 
cavity  surface.  We  then  show  that  two-dimensional  gaseous  diffusion 
models  cannot  always  be  used  to  predict  gaseous  diffusion  across 
axisymmetric  cavities.  This  naturally  led  us  to  the  formulation  of 
convective  diffusion  of  gases  for  axisymmetric  cavity  shapes.  We 
solve  the  steady  diffusion  equation  in  cylindrical  coordinates  as  it 
applies  to  the  gaseous  diffusion  across  an  axisymmetric  cavity  surface. 
Although  the  general  solution  of  the  governing  differential  equation 
is  fairly  straightforward,  the  peculiar  set  of  boundary  conditions 
for  the  present  idealized  representation  of  the  flow  and  cavity  of 
Interest  here  makes  the  problem  challenging. 

The  solution  of  the  problem  is  built  up  by  superposition  starting 
from  the  known  point  source  solution  in  a  moving  fluid.  First,  we  use 
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the  point  source  solution  to  derive  the  solution  to  the  governing  ] 

differential  equation  for  a  “unit  ring  source."  Second,  we  integrate  I 

the  ring  source  of  unknown  strength  over  the  length  of  the  doubly 
infinite  circular  cylinder  in  the  strearawise  direction.  The 
application  of  boundary  conditions  on  this  cylindrical  cavity  surface 
results  in  a  system  of  integral  equations  which  must  be  solved  in  order 
to  determine  the  source  strength.  From  this  solution,  the  mass  flow 
rate  of  gases  diffusing  across  the  cavity  may  be  obtained. 

As  pointed  out  by  Greitzer  (20),  there  is  a  general  lack  of 
theoretical  and  experimental  results  available  in  the  literature  that 
characterize  the  unsteady  cavity  as  an  active  dynamic  element. 

Therefore,  an  effort  was  made  to  understand  the  behavior  of  unsteady 
cavities  experimentally  in  this  study.  We  conducted  experiments  on  a 
cavitating  wedge  foil  oscillating  in  pitch.  Specifically,  we  were 
interested  in  measuring  the  cavity  lengths  and  cavity  pressures  (steady 
and  nonsteady)  as  a  function  of  oscillation  frequency  and 
characterizing  the  source-like  nature  of  oscillating  cavities.  We 
found  that  an  oscillating  cavity  has  the  behavior  of  a  dipole  source 
and  that  the  unsteady  cavity  pressures  are  significant  at  low 
frequencies  of  oscillation.  We  also  found  that  the  motion  of  the 
cavity  is  far  from  being  sinusoidal  even  though  the  motion  of  the  foil 
was  sinusoidal.  In  particular,  the  first  harmonic  overtone  was  found 
to  be  significant.  In  Chapter  3  the  details  of  this  experimental 
Investigation  are  presented  and  the  results  are  discussed. 


8 


The  experimental  results  presented  in  Chapter  3  as  well  as  those 
of  Jiang  and  Leehey  (17)  have  indicated  clearly  the  effect  of  non¬ 
condensable  gas  on  the  behavior  of  an  unsteady  cavity.  While  several 
investigations  in  the  past  have  addressed  the  effect  of  steady  gaseous 
diffusion  into  the  cavity,  there  is  very  little  literature,  if  any, 
that  deals  with  the  analysis  of  unsteady  diffusion  across  cavity 
surfaces.  In  Chapter  4,  we  formulate  and  solve  the  two-dimensional 
unsteady  diffusion  problem.  The  problem,  as  formulated,  accounts 
for  both  the  cavity  pressure  fluctuations  and  cavity  length 
fluctuations.  The  solution  of  the  unsteady  diffusion  equation  is 
obtained  by  the  separation  of  variables.  The  expression  for  mass 
diffusion  rate  is  obtained  in  a  manner  similar  to  the  method  outlined 
by  Parkin  (8)  for  steady  two-dimensional  flow.  The  result  for 
unsteady  mass  diffusion  rate  across  a  two-dimensional  cavity  surface 
matches  the  steady  state  solution  predicted  by  Parkin  (8)  in  the 
limit  when  the  frequency  of  cavity  oscillation  vanishes. 

In  Chapter  5,  we  consider  the  two  dimensional  steady  cavity 
closure  as  well  as  nonsteady  cavity  closure  condition  when  the  cavity 
is  not  purely  vaporous.  The  steady  closure  condition  is  seen  to  be 
rather  straightforward  to  formulate  and  analyze.  In  analyzing  the 
nonsteady  gaseous  cavity  closure,  we  find  that  it  is  important  to 
consider  the  phase  differences  that  exist  between  the  body  motion  and 
i)  cavity  volume  (or  almost  equivalently  the  cavity  length) 
fluctuation,  ii)  the  unsteady  gaseous  diffusion  through  the  entire 
cavity  surface,  and  iii)  gas  pressure  fluctuation  within  the  cavity 
respectively.  In  the  absence  of  any  available  theoretical  or 
experimental  data  on  the  phase  lag  between  the  body  motion  and  cavity 
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length  motion,  we  appeal  to  a  physically  driven  argument  that  cavity 
surface  waves  originating  from  the  upper  surface  and  lower  surface  of  a 
foil  must  terminate  at  the  same  point  at  the  same  instant.  This 
results  in  the  phase  lag  being  numerically  equal  to  the  modified 
reduced  frequency  of  foil  oscillation.  The  nonsteady  gaseous  cavity 
closure  is  then  formulated  using  the  nonsteady  continuity  equation. 

With  the  help  of  the  experimental  data  reported  above,  continuity 
enables  the  determination  of  magnitude  and  phase  of  the  nonsteady 
entrained  gas  behind  the  cavity  in  terms  of  the  reduced  frequency  of 
oscillation  k'f  the  nonsteady  mass  diffusion  through  the  cavity 
surface,  the  cavity  pressure  fluctuations  and  the  cavity  volume 
fluctuations. 

The  reader  will  find  that  an  analytical  approach  has  been  used 
throughout  the  present  study.  The  reason  for  this  being,  analytical 
solutions  render  better  physical  insight  than  numerical  solutions. 
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CHAPTER  2 

GASEOUS  DIFFUSION  STUDIES  ACROSS  AN  AXI SYMMETRIC 
PERMEABLE  CAVITY  SURFACE 


2. 1  Introduction 

A  thorough  analysis  of  the  gaseous  diffusion  process  across  a 
cavity  surface  is  fundamental  to  the  understanding  of  both  natural  and 
ventilated  cavity  flows  because  diffusion  into  or  out  of  the  liquid 
can  influence  the  cavity  pressure  to  a  great  extent.  Analysis  of 
gaseous  diffusion  across  axisymmetric  cavity  surfaces  is  of  particular 
Interest  because  of  the  frequent  occurrence  of  such  cavity  flows  in 
marine  applications. 

Steady  diffusion  of  gases  across  a  two-dimensional  cavity  has 
been  studied  by  Parkin  (8).  Parkin  used  this  planar  solution  to 
approximate  the  rate  of  molecular  diffusion  into  a  microbubble. 

Brennen  (9)  considered  turbulent  gaseous  diffusion  into  a  large 
two-dimensional  cavity  by  using  a  model  which  assumes  that  the 
turbulence  is  associated  with  the  turbulent  shear  layer  on  the  cavity 
surface.  Billet  and  Weir  (13)  have  attempted  to  use  Brennen's  model 
to  predict  the  rate  of  diffusion  across  axisymmetric  cavity  surfaces. 
They  found  when  Brennen’s  two-dimensional  solution  is  extended  to  the 
axisymmetric  case  by  wrapping  the  two-dimensional  solution  around  the 
cavity,  that  the  observed  diffusion  rate  exceeds  the  analytical 
prediction  by  nearly  one  order  of  magnitude.  They  circumvented  this 
discrepancy  by  using  their  experimental  data  to  choose  a  constant 
multiplicative  factor  in  Brennen's  model. 
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Brennen's  solution  assumes  that  the  diffusion  occurs  through  the 
turbulent  eddies  within  the  turbulent  boundary  layer  on  the  cavity 
surface.  Hence  Brennen's  solution  is  applicable  only  where  the 
boundary  layer  has  had  a  chance  to  grow  on  the  surface  of  the  body 
before  a  cavity  is  formed,  as  for  example,  behind  the  trailing  edge 
of  a  hydrofoil.  But  in  cases  where  the  cavity  emerges  from  the 
leading  edge  of  a  hydrofoil  at  small  angle  of  attack  or  the  cavity  is 
on  a  zero  caliber  ogive  as  in  the  experiments  of  Billet  and  Weir, 
the  boundary  layer  on  the  cavity  near  the  leading  edge  is  laminar 
and  very  thin;  and  since  the  cavity  surface  is  then  a  laminar  free 
shear  layer,  the  rate  of  growth  of  this  part  of  the  boundary  layer 
is  small.  For  such  flows  it  appears  as  though  a  different  approach 
has  to  be  followed  in  modeling  the  turbulent  diffusion  phenomenon. 

In  many  cases,  this  laminar  shear  layer  is  of  short  axial  extent 
and  the  shear  layer  is  turbulent  thereafter. 

In  this  chapter,  an  attempt  is  made  to  clear  some  of  the 
discrepancies  observed  by  Billet  and  Weir  (13)  in  their  effort  to 
compare  their  measurements  with  Brennen's  theory.  A  convective  slug 
flow  model  is  first  used  to  estimate  the  maximum  possible  mass 
convective  diffusion  rate  of  gas  through  the  turbulent  boundary 
layer  on  the  axisymmetric  cavity  surface.  These  estimates  are  then 
compared  with  the  experimental  results  of  Billet  and  Weir  (13)  and 
Brennen's  theory.  Then,  a  turbulent  diffusivity  (that  is 
representative  of  free  shear  flows)  is  used  instead  of  molecular 
diffusivity  as  in  Parkin's  model,  to  compare  with  the  experimental 
results.  A  mathematical  model  is  then  formulated  for  gaseous 
diffusion  across  an  axisymmetric  cavity  surface. 
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2.2  Convective  Slug  Flow  Model  for  Mass  Diffusion  Rate  Estimation 
Brennen's  theory  assumes  that  all  of  the  diffusion  occurs  from 
turbulent  eddies  that  are  within  the  boundary  layer.  For  cavity 
flows  wherein  the  cavity  emerges  from  the  leading  edge,  the  boundary 
layer  on  the  cavity  surface  is  very  thin  and  Brennen's  theory  may 
not  quite  be  applicable.  In  order  to  demonstrate  whether  or  not 
Brennen's  theory  is  applicable  in  such  situations  and  to  assist 
speculations  about  other  possible  mechanisms  for  gaseous  diffusion 
in  cavity  flows,  the  following  conservative  estimates  are  carried 
out  using  a  slug  flow  model. 

We  assume  that  the  quantity  of  dissolved  gases  available  for 
diffusion  is  contained  within  a  volume  of  liquid  determined  by  the 
maximum  thickness  of  the  turbulent  boundary  layer  and  the  length  of 
the  cavity  as  shown  in  Figure  2.1.  We  also  assume  that  all  gas  in 
this  volume  diffuses  into  the  cavity  instantaneously  and  then  the 
slug  of  volume  in  the  shape  of  a  hollow  cylinder  of  mean  radius 
"a,”  length  "L"  and  thickness  "5,"  equal  to  the  maximum  shear  layer 
thickness  at  the  cavity  terminus,  is  swept  downstream  and  a  fresh 
supply  of  gases  becomes  available  every  L/U  seconds. 

If  we  denote  by  c,  the  concentration  difference  across  the 
cavity  surface,  the  maximum  possible  diffusion  rate  is  estimated 
as 


or 


M  =■  (2;raL5)c(U  /L) 


M  *  2waU  6c 


(2.1) 
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TURBULENT  BOUNDARY  LAYER 


Figure  2.1.  Slug  Flow  Model  for  Estinating  the  Gaseous 
Diffusion  across  an  Axisysretric  Cavity. 


The  mass  diffusion  rates  predicted  by  the  slug-flow  model  are  shown 


in  Table  2.1  for  cavities  behind  zero  caliber  ogive  models.  Also 


shown  in  the  table  are  the  experimental  values  of  mass  diffusion 


rates  from  Billet  and  Weir  (13)  and  the  mass  diffusion  rates  predicted 


by  Brennen's  theory  (9).  A  comparison  of  the  diffusion  rates 


predicted  by  the  slug-flow  model  and  the  experiments  shows  that  even 
by  conservative  estimates,  the  maximum  gas  diffusion  possible  through 


the  turbulent  boundary  layer  is  nearly  one  order  of  magnitude  less 


than  the  experimentally  observed  values, 


2.3  Modified  Parkin  Model 


Following  the  analysis  in  Ref.  (8),  the  mass  diffusion  rate 


across  a  two-dimensional  cavity  surface,  with  the  molecular 


diffusivity  <  being  replaced  by  a  turbulent  diffusivity  vt  is  given 


4Pvtc  r  _  /-  , 

-  [/2U  -  ]  , 

✓it  1 


(2.2) 


where 


mass  flow  rate  per  unit  cavity  width 


density  of  water 


turbulent  diffusivity 


concentration  difference 


U«a/2vt 


length  of  two-dimensional  cavity. 


Wrapping  the  two-dimensional  solution  around  an  axisymmetric  cavity 
of  radius  "a,"  one  can  write  the  mass  flow  rate  across  an  axisymmetric 
cavity  as 


M  »  8pvtc  /ir  a  [/2A£  -  ] 


(2.3) 


For  a  zero  caliber  ogive,  the  non-dimensional  cavity  length  £  and 
cavity  radius  is  related  to  the  cavitation  number  a  by  the  relations* 


£ 


and 


a 


0.715 

0.34 

a 


D  . 


(2.4a) 


(2.4b) 


In  Eq.  (2.4b),  D  represents  the  diameter  of  the  axisymmetric  ogive 
body.  With  Eqs.  (2.4)  substituted  in  Eq.  (2.3), 

M  -  10.14pcDvto0,34  [l.382A1/2<f 1/2  -  0.886]  .  (2.5) 


We  now  devise  a  turbulent  diffusivity  parameter  vt  characteristic 
of  a  free  shear  layer  as 


Vt  ~ 


(2.6) 


where 


£ 

m 


*  r.ra.s.  turbulence  level 
■  mixing  length. 


*(See,  for  example,  Billet,  M.  L. ,  J.  W.  Holl  and  D.  S.  Weir, 
"Geometric  Description  of  Developed  Cavities  on  Zero  and 
Quarter  Caliber  Ogive  Bodies,"  ARL/PSU  TM  74-136  [1974].) 
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For  a  planar  mixing  layer* 

£  ~  0.075 
01 

and 


~  0.09U  , 

OB  9 


where  5  is  the  layer  thickness.  Therefore  we  will  take 


(2.7a) 


(2.7b) 


vt  -  0.0063  6U.y  .  (2.7c)  ] 

1 

In  the  following  comparison,  the  constant  of  proportionality  y  is 
taken  to  be  unity.  Figures  2.2  and  2.3  show  the  variation  of 

i 

volume  flow  rates  predicted  by  Eq.  (2.5)  and  the  corresponding  j 

experimental  values  of  Billet  and  Weir  (13).  It  is  observed  that  j 

the  theoretical  and  experimental  values  compare  well;  at  least  they 

are  of  the  same  order  of  magnitude!  This  good  agreement  is  to  be 

expected  if  it  is  recognized  that  near  the  leading  edge  (where  the  ' 

I 

concentration  gradients  are  very  high)  the  boundary  layer  thickness  j 

I 

is  very  small  and  thus  extension  of  a  two-dimensional  result  to 
axisymmetric  cases  appears  to  give  a  good  approximation  even  though 
the  concentration  layer  is  comparable  to  the  cavity  radius  at  the 
tail  end  of  the  cavity  (see  Figure  2.4  below). 


*See,  for  example.  Launder,  B.  E.  and  D.  B.  Spalding,  Mathematical 
Models  for  Turbulence,  Academic  Press  [1972). 
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Flaure  2.2  Flow  Coefficient  versus  Cavitation  Number,  -  1830  cm/sec. 
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Coefficient  versus  Cavitation  Number,  U  ”  910  cm/sec 
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2.4  Mathematical  Analysis 

The  discrepancy  that  exists  between  the  experimental  results 
and  existing  mathematical  models  for  gaseous  diffusion  across 
axisymmetric  cavities  prompted  us  to  look  at  the  validity  of 
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extending  existing  two-dimensional  diffusion  models  to  axisymmetric 
diffusion  studies.  While  it  is  true  that  two-dimensional  solutions 
can  be  extended  to  include  axisymmetric  cases,  such  extensions  hold 
good  only  when  the  concentration  boundary  layer  thickness  on  the 
cavity  is  small  compared  to  the  radius  of  cavity.  This  is  a  rather 
serious  restriction  for  flows  normally  encountered  where  the 
concentration  layer  thickness  is  not  necessarily  small,  relative  to 
the  radius  of  the  cavity.  Just  to  exemplify  this  point,  the 
concentration  layer  profile  at  the  end  of  the  cavity  has  been  plotted 
in  Figure  2.4,  the  flow  conditions  being  representative  of  experi¬ 
ments  reported  in  Ref.  (13).  It  is  clear  that  the  concentration 
layer  thickness  is  comparable  to  the  radius  of  cavity.  Thus,  a 
more  accurate  mathematical  representation  for  gaseous  diffusion 
across  an  axisymmetric  cavity  surface  is  clearly  warranted. 

Gaseous  diffusion  will  occur  when  there  is  a  dissolved  gas 
concentration  gradient  between  the  free  stream  and  the  liquid  on  the 


cavity  surface.  If  is  the  measured  dissolved  gas  content  in 


parts  per  million  (ppm)  by  moles,  then  by  Henry's  law,  the  saturation 
partial  pressure  of  the  gas  in  the  free  stream  is  given  by 


PFS  "  V 


(2.8) 


j; 
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where  8  is  the  Henry's  law  constant  which  depends  on  temperature.  If 
p  is  the  partial  pressure  of  non-condensable  gas  in  the  cavity, 
then  the  mean  concentration  difference  expressed  in  ppm  by  moles  is 


c  "  °i  "  T 


(2.9) 


Consider  an  axisymmetric  flow  of  water  containing  dissolved  gases 
which  has  the  concentration  C(X,R)  at  any  point  in  the  region  R  >  a, 
where  "a”  is  the  radius  of  the  axisymmetric  cavity.  The  undisturbed 
flow  is  taken  to  have  a  constant  velocity  U  in  the  positive 

OB 

X-direction.  In  the  interest  of  simplicity,  all  momentum  boundary 
layer  perturbation  velocities  u  and  v  in  the  X  and  R  directions 
respectively  are  neglected.  The  free  stream  velocity  applies  to 
the  entire  flow  field.  The  mass  diffusion  of  gases  through  the 
permeable  axisymmetric  cavity  surface  occurs  between  X  -  0  and  X  *  L 
as  shown  in  Figure  2.5.  The  remaining  portion  of  the  X-axis  is 
supposed  to  be  impermeable  to  gas  diffusion.  On  the  permeable 
surface,  let  the  concentration  have  the  constant  value  and  the 
concentration  of  dissolved  gases  at  points  far  from  the  permeable 
surface  have  the  value  C^.  Then,  the  concentration  difference  c(X,R) 
is  defined  as 


c(X,R)  -  C(X,R)  -  C 


(2.10) 
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Let  tc  be  the  mass  diffusivity  parameter.  We  do  not  identify  the 


nature  of  the  mass  diffusivity  parameter  as  either  molecular  or 


turbulent.  In  either  case,  the  governing  differential  equation  for 


the  concentration  difference  has  the  form  [see,  for  example, 


Carslaw  and  Jaeger  (22)3 


t,  2 

U  —  =■  icV  c  , 
®  3x  ’ 


(2.11) 


where  *7  is  the  Laplacian  operator.  By  non-ditnensionalizing  the 


linear  distances  along  X,  Y,  Z  and  R  by  the  radius  of  the  cavity 


such  that  x  *  X/a,  y  ■  Y/a,  z  =  Z/a  and  r  =■  R/a,  we  can  write 


Gq.  (2.11)  as 


3c  _  k  _2 
3x  aU  C 


(2.12) 


In  order  to  determine  the  solution  of  Eq.  (2.12)  subject  to  the 


prescribed  boundary  conditions,  we  first  write  the  solution  of 


Eq.  (2.12)  when  a  source  of  strength  ”q”  is  located  at  the  origin 


in  a  moving  fluid  (22),  as 


c(x,y,z) 


-A(r ’  -  x) 


4irkr  a 


(2.13) 


where 


q  ■  mass  flux  of  gas 


k  ■  mass  conductivity  *  icp 


r'  *  Vx^  +  y^  +  z^  ,  the  distance  between  the 


field  point  and  the  origin 


k  ■  mass  diffusivity 


p  *  density  of  the  liquid 


\  *  Ua/2ic  ,  a  dimensionless  parameter. 
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If  the  source  is  located  at  (5,n,C)  instead  of  at  the  origin, 
the  solution  of  Eq.  (2.12)  will  be  modified  as  follows: 

,  _  .  q  -X(r ' '  -  x)  , 

c(x,y,z;  5,n,s)  =— — n~  e  »  (2.14) 

4irkr  a 

1 1  j  2  2  2 

where  r  =■  V(x  -  ?)  +  (y  -  n)  +  (z  -  5)  ,  is  the  distance  between 

the  source  point  (5,11,5)  and  the  field  point  (x,y,z). 

We  now  consider  a  "ring  source"  of  unit  radius  and  unit  strength 
located  at  x  «  The  plane  of  the  ring  is  perpendicular  to  the 
x-axis  and  the  x-axis  passes  through  the  center  of  the  ring.  Consider 
an  element  of  the  ring  d0,  as  shown  in  Figure  2.6.  Then  the  strength 
of  this  element  will  be  d0/2ir  and  from  Eq .  (2.14),  the  concentration 
at  a  point  (x,r)  due  to  the  flux  of  concentration  at  (5,1,0)  will  be 


6c(x,r ; 5, 1 , 0)  - 


d0  -X(r  -  (x  -  5)) 

q  2.  1 »  6 

8ir  kr  a 


(2.15) 


Integrating  the  right-hand  side  of  the  above  equation  between  the 
limits  of  0  and  2ir,  we  obtain  the  concentration  at  a  field  point 
(x,r)  due  to  a  unit  ring  source  at  x  =  5  as 

X(x  -  5)  2ir  -Xr 

c(x,r;  5)  =  “ - 5 -  /  ~r ,  ,  •  d0  ,  (2.16 

8ir  ka  0 

where  the  distance  r''  is  now  given  by 

r ' '  =  V(x  -  +  r^  +  1  -  2r  cos  0  (2.17 

The  integral  in  Eq.  (2.16)  has  been  evaluated  in  Appendix  A.l  for 
large  values  of  X  (of  the  order  2000)  by  using  Laplace's  method  and 


the  result  for  the  concentration  is 


c(x,r ; £)  =< 


eX(x-£)  e-X  Y(x  -  £)  +  (r  -  1)‘ 

8ir3/2pa3/2/CTir  [(x  -  £) 2  +  (r  -  l)2] 


2 1 1/ 4 


.  (2.1! 


Equation  (2.18)  has  the  desired  properties  relevant  to  the  physical 
problem,  viz.,  c  ♦  0  as  x  +  -  •  and  c  +  0  as  r  ♦  «. 

Now  suppose  that  the  mass  flux  per  unit  length  between  x  =  -  «• 
and  x  *  +  «•  is  given  by  the  function  f(x).  Then  the  function  c(x,r) 
can  be  expressed  in  terms  of  this  source  strength  per  unit  length  as 

A  “  p_x  V(x  -  £ )2  +  (r  -  l)2 

c(x,r)  ■  ~  /  f(£)eX  5)  [ - r -  ]d£  ,  (2.1 

✓r  —  ((x  -  O  +  (r  -  l)zJl/ 


where 


-3/2  3/2,— 

8it  pa  yk  U 


In  the  integral  in  Eq.  (2.19),  the  limit  -  «•  to  +  »  appears  instead 
of  0  to  l  because  the  unit  ring  source  solution  does  not  satisfy  the 
boundary  condition  at  r  *  1:  dc/dr|  »  0.  If  we  now  require  that 
c(x,l)  ®  cn  *  Cn  -  C  ,  in  the  interval  0  <  x  <  l,  and  dc/dr|  ,  »  0 
in  the  intervals  -  «  <  x  <  0  and  4  <  x  <  »,  the  source  strength  f(x) 
is  determined  by  the  three  integral  equations  of  the  form 

x  -2X(£  -  x  ) 

c0  0  f,(£)  2  f -(£)  *  f-(£)e 

r-/  -J — *  f  -j — dg +  i  -2 — - — d5 

_oo  /x^  -  £  0  /x^  -  £  x2  /£  -  x^ 


/ £  -  x. 


-2X(£  -  xj 


It,' 


•WwS 


« 

k 


2X(5  -  xL) 


-2A U  ~  x, ) 

1  f„<5)e  1 


2A(5  -  x L) 

- d?  , 

;  "X1 

1  1  f,U) 

X3  f,(o 

where  the  positive  branch  of  the  square  root  has  been  taken.  The 
functions  f^,  f^  and  f^  are  the  source  strength  distributions  in  the 
intervals  -  <■  <  <  0,  0  <  ^  <  l  and  l  <  <  «  respectively.  The 

three  source  functions  are  taken  to  be  zero  outside  their  specified 
intervals . 

The  three  integral  equations  (2.21)  to  (2.23)  are  coupled  and 
rather  difficult  to  solve  in  their  present  form.  In  order  to  make 
the  solution  tractable,  we  simplify  some  of  the  integrals  appearing 
in  the  three  equations  (see  Appendix  A. 2)  recognizing  the  fact  that 
the  nondiraensional  number  A  is  fairly  large.  Thus,  the  simplified 
integral  equations  are 


J  —=<n  +  f  7==d*+  ixW  * 


(2.21a) 


^  f.<5> 

r  1 


/x^  -  5 


d5+ 


(2.22a) 


0  f,(5)  f,(5)  3  f,(C) 

/ - -1 - dg  +  f  -  dC  »  f  -  ■—  dg  +  f.(x  )  fr-  -  0  . 

-co  /x-j  -  5  o  ^x3  ~  5  £  ^x3  -  € 


(2.23a) 

The  integral  equation  (2.22a)  has  been  solved  using  Picard's  method* 
in  Appendix  A3  and  the  result  for  f  is 


ZAX. 

fjxj  -ct.  *,  <  0  , 


(2.24) 


and  c^  is  an  arbitrary  constant.  The  function  f^(x^)  may  now  be 
substituted  in  the  integral  equation  (2.21a)  and  the  function 
f,(x  )  obtained  by  the  Laplace  transform  method  (see  Appendix  A. 4) 


f2(x2)  '  f^+cl)e  2erfc/21x:  -  -  1  -  ,  (2.25] 

/2X1TX2 


where  the  parameter  X  has  been  defined  as  part  of  Eq.  (2.13),  the 
quantity  A  is  given  by  Eq.  (2.20)  and 


*See,  for  example,  Hildebrand,  Methods  in  Applied  Mathematics, 
Prentice  Hall  Inc.  N.  Y.  1952. 
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where  erfc(t)  denotes  the  complimentary  error  function 


2  r  -x 

erfc(t)  »  —  /  e  dx  .  (2 

/ir  c 

Having  found  f^(x^)  and  £2^2) »  one  can  find  the  function  f^x^) 
from  Eq.  (2.23a).  We  have  not  given  this  last  solution  here  even 
though  the  way  to  find  it  seems  clear.  The  reason  is  that  £2^2) 
controls  the  diffusion  into  the  cavity  and  f^  is  not  required  to 
determine  f2*  We  note  also  that  the  source  strength  is  strongest 
near  the  cavity  separation  point  so  most  of  the  mass  flux  will 
occur  at  the  start  of  the  cavity  and  not  near  the  downstream  end. 

The  rate  at  which  mass  is  diffused  along  the  entire  length  of 
the  axisymmetric  cavity  is 


*  dM  r  *  r  'u 

"  ■  jr  -  /  • 

0 


fn  . 


^erfc/2 \l  .  12 1  1 


ft,1 


e2A*erfc/2U  ~  —  f— — - 1-  +  ...]  ,  (2.29) 
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in  which  case, 
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(2.30) 


If  we  assume  that  the  constant  c^  has  the  same  order  of  magnitude 
2 

as  8irpa  U^Cq,  w®  can  conclude  from  Eq.  (2.30)  that  the  expression 
for  mass  flow  rate  across  an  axisymmetric  cavity  is  the  same  as 
the  two-dimensional  wrap  around  solution.  This  conclusion,  of 
course,  is  based  on  the  simplifications  we  have  effected  in  the 
three  coupled  integral  equations  (2.21),  (2.22)  and  (2.23). 

These  same  approximations  were  employed  by  Parkin  in  his  theory. 
Nonetheless,  the  good  correlations  observed  between  Billet-Weir 
experimental  results  and  modified  Parkin's  model  (see  Figures  2.2 
and  2.3)  support  the  conclusions  that  we  draw  from  the  axisymmetric 
cavity  diffusion  analysis.  The  reason  for  this  being  that  most  of 
the  diffusion  takes  place  across  the  upstream  parts  of  the  cavity 
where  the  diffusion  layer  is  still  thin  and  wrap-around  is  still 


useful 


2.5  Recapitulation 

We  have  explored  some  of  the  aspects  of  gaseous  diffusion  into 
axisymmetric  cavities.  A  modifed  Parkin  model  may  be  better  suited 
for  the  prediction  of  gaseous  diffusion  across  cavities  that  emerge 
from  the  leading  edge  than  Brennen’s  model.  We  have  demonstrated 
through  an  axisymmetric  gaseous  diffusion  analysis  that  (within  the 
assumptions  and  approximations)  two-dimensional  gaseous  diffusion 
solutions  can  be  extended  to  the  axisymmetric  cavity  with  the 
result  that  the  turbulent  mixing  length  diffusivity  of  Launder  and 
Spalding  will  be  useful  for  the  determination  of  diffusive  mass 
flow  into  the  two-dimensional  flows  studied  below. 
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CHAPTER  3 


EXPERIMENTAL  INVESTIGATION  OF  THE  BEHAVIOR  OF  NONSTEADY 
CAVITIES  BEHIND  A  WEDGE  OSCILLATING  IN  PITCH 


3.1  Introduction 

An  unsteady  cavity  flow  about  a  solid  body  in  a  pure  liquid  may 
be  represented  mathematically  with  the  help  of  cavity  flow  theory  and 
it  can  also  be  understood  physically.  The  presence  of  dissolved  gases 
in  the  liquid  makes  the  same  phenomenon  difficult  to  understand.  A 
cavity  in  a  pure  liquid  is  filled  only  with  the  liquid  vapor.  So  the 
volume  of  an  unsteady  cavity  filled  only  with  vapor  grows  and 
collapses  in  volume  according  to  the  motion  of  the  body  or  according 
to  the  imposed  pressure  field.  In  the  absence  of  other  flow  distur¬ 
bances,  the  vapor  is  absorbed  and  released  almost  instantaneously 
because  the  time  scale  of  the  process  of  vapor  release  and  absorption 
is  extremely  small  compared  to  the  time  scale  of  cavity  motion  and 
hence  cavity  vapor  pressure  fluctuations  are  not  seen.  In  such  a 
case,  the  pressure  inside  the  cavity  equals  that  of  vapor  pressure 
of  the  ambient  liquid  at  the  equilibrium  temperature  at  all  times, 
provided  the  temperature  is  not  too  near  the  critical  temperature. 

This  situation,  however,  does  not  occur  in  the  laboratory  or 
in  practice.  There  are  always  certain  amounts  of  free  and  dissolved 
gases  present  in  the  liquid.  On  a  cavitating  body,  these  gases 
diffuse  across  the  cavity  surface  and  are  entrained  by  the  turbulent 
flow  at  the  trailing  edge  of  the  cavity.  This  diffusion-entrainment 
process  is  thought  to  influence  cavity  volume  fluctuations  and 
associated  cavity  pressure  fluctuations  in  an  unsteady  cavity. 


Several  investigations  have  been  devoted  in  the  past  to  the 


analysis  of  steady  gaseous  diffusion  across  the  cavity  surface. 

Parkin  and  Kermeen  (8)  studied  the  molecular  gaseous  diffusion  into 
a  microbubble,  whereas  Brennen  (9)  considered  gaseous  diffusion  from 
a  turbulent  flow  into  a  large  cavity.  Very  few  unsteady  cavity  flow 
studies  (both  analytical  and  experimental)  have  been  done  in  the  past 
in  which  cavity  pressure  fluctuations  are  considered,  largely  due  to 
the  complexities  involved  in  such  studies.  Although  Leehey  and 
Jiang  (16)  measured  the  unsteady  cavity  pressures  behind  a  three- 
dimensional  elliptical  foil  oscillating  in  pitch,  there  appears  to 
be  no  experimental  study  available  aimed  at  a  better  understanding 
of  the  behavior  of  unsteady  cavities,  particularly  (a)  the  cavity 
pressure  fluctuations  and  their  dependence  on  the  frequency  of 
oscillations  of  the  body;  (b)  the  source-like  nature  of  the  unsteady 
cavity;  and  (c)  the  cavity  length  fluctuation  as  a  function  of 
frequency.  The  present  experimental  study  is  aimed  at  understanding 
these  aspects  of  unsteady  two-dimensional  cavities. 

3.2  Experimental  Apparatus 
3.2.1  Six-Inch  Water  Tunnel 

The  experiment  was  conducted  in  the  6-inch  recirculating 
variable  pressure  water  tunnel  at  ARL/PSU  (21).  In  this  tunnel,  the 
test  section  velocity  was  variable  between  0  and  50  feet  per  second 
and  the  static  pressure  was  variable  continuously  between  0.5  psia 
and  50  psia.  The  velocity  was  indicated  on  a  digital  voltmeter 
which  read  the  pressure  drop  between  the  settling  chamber  and  the 


test  section.  The  test  section  static  pressure  was  read  by  a  Heise 
type  Bourdon-tube  pressure  gage.  The  test  section  was  24  inches 
long  and  was  made  of  lucite  to  enable  visual  observations  and  video 
recordings  to  be  made  of  the  unsteady  cavity  flow  phenomena. 

3.2.2  Model  and  Shaft 

The  foil  used  in  this  test  was  a  wedge  of  rectangular  planform 
measuring  two  inches  wide  and  four  inches  long  and  was  made  of 
stainless  steel.  A  small  slot  was  made  on  the  cavitating  surface  of 
the  foil  to  mount  a  pressure  transducer  that  would  measure  the 
unsteady  cavity  pressures.  A  small  hole  was  made  in  the  shaft  running 
axially  to  the  middle  of  the  foil  to  lead  the  pressure  transducer 
wires  and  also  a  static  cavity  pressure  tap  out  of  the  tunnel.  The 
support  shaft  was  made  of  stainless  steel  and  the  middle  portion  of 
the  shaft  formed  the  trailing  section  of  the  wedge  foil.  The  foil  and 
shaft  were  fastened  together  by  a  pair  of  dowel  pins  and  set  screws. 
Two  views  of  the  foil-shaft  assembly  mounted  in  the  test  section  along 
with  the  hydraulic  drive  are  shown  in  Figure  3.1. 

3.2.3  Drive  System 

The  drive  system  used  was  essentially  a  modified  control  surface 
actuator.  It  consisted  of  two  plungers  actuated  by  a  hydraulic 
servo.  The  plungers  were  connected  to  the  pitch  drive  plate  fastened 
to  the  shaft  by  means  of  a  wire  rope.  When  the  servo  was  actuated 
by  an  amplified  sine  wave  signal,  the  shaft  executed  nearly  sinusoidal 
motion  in  pitch.  Because  of  the  slider  crank  mechanism,  though  the 


Foil  Shaft  Assembly . 
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motion  of  the  plunger  was  sinusoidal,  the  motion  of  the  shaft  was  not. 
However,  an  analysis  of  the  error  signals  from  the  feed  back  unit 
showed  higher  harmonics  whose  amplitudes  were  at  least  40  dB  lower 
than  the  fundamental  frequency  of  oscillation.  The  angle  of  attack 
of  the  foil  was  varied  by  changing  the  d.c.  offset  to  the  hydraulic 
servo.  The  frequency  response  of  the  pitch  drive  assembly  which  is 
flat  up  to  about  10  Hertz  is  shown  in  Figure  3.2.  The  static  angle 
of  attack  of  the  foil  was  variable  continuously  from  zero  to  twenty 
degrees.  The  dynamic  angle  of  attack  amplitude  was  variable 
continuously  from  zero  to  fifteen  degrees  peak  to  peak.  Using  this 
setup,  the  static  angle  of  attack  could  be  adjusted  to  an  accuracy 
of  0.02  degrees  and  the  amplitude  of  foil  oscillation  could  be 
adjusted  to  an  accuracy  of  ±  0.05  degrees.  Figure  3.3  shows  a 
typical  cavity  behind  the  wedge  foil. 

3.2.4  Instrumentation 

Unsteady  cavity  pressure  was  measured  by  mounting  a  crystal 
transducer  flush  with  the  cavitating  surface  of  the  foil.  This  Barium 
Titanate  crystal  transducer  was  calibrated  in  a  slosh  tube  (22)  with 
a  LC-10  hydrophone  as  the  reference.  The  frequency  response  of  the 
crystal  transducer  mounted  on  the  foil  is  shown  in  Figure  3.4.  The 
static  pressure  in  the  cavity  was  measured  by  a  pressure  tap  leading 
to  the  cavity  from  a  Validyne  diaphragm  type  pressure  transducer. 

The  unsteady  pressures  upstream  and  downstream  of  the  foil  were  also 
measured  with  a  Validyne  pressure  transducer.  This  transducer  was 
suspended  from  an  independent  support  to  eliminate  possible  effects 
of  tunnel  vibration  on  the  pressure  transducer  output. 


The  amount  of  free  air  in  the  tunnel  water  was  determined  by 
measuring  the  speed  of  sound  in  the  upstream  part  of  the  test  section. 
This  was  achieved  by  transmitting  a  gated  high  frequency  pulse  (about 
120  kHz)  through  a  LC-10  hydrophone  and  receiving  it  diametrically 
across  the  test  section  by  a  crystal  transducer  (Barium  Titanate 
crystal)  as  shown  in  Figure  3.5.  The  distance  between  the  two  trans¬ 
ducers  for  the  sound  speed  measurement  was  7.443  inches.  By  measuring 
the  time  delay  between  the  transmitted  and  received  pulse  on  a  high 
speed  oscilloscope  (Tektronix  7633),  the  speed  of  sound  was 
determined.  The  accuracy  of  speed  measurement  is  within  i  30  ft/sec. 
An  example  of  the  sound  speed  calculation  and  the  estimation  of  free 
air  in  the  tunnel  water  is  given  in  Appendix  B.2.  The  dissolved  air 
in  the  tunnel  water  was  measured  by  the  conventional  Van  Slyke 
Apparatus.  During  the  course  of  the  experiment,  the  dissolved  air 
content  ranged  between  7  ppm  and  16  ppm  by  moles  while  the  free  air 
ranged  from  3  ppm  (speed  of  sound  =  4490  ft/sec)  to  10  ppm  by 
volume  (speed  of  sound  =  3730  ft/sec).  The  measurement  of  free  gas 
and  dissolved  gas  in  the  tunnel  ensured  that  there  was  enough  gas 
supply  in  the  tunnel  water  to  maintain  dissolved  gas  diffusion  and 
entrainment  at  all  times . 

All  the  unsteady  pressure  data  were  analyzed  through  a  real  time 
spectral  analyzer  (Spectral  Dynamics  SD-300  and  SD-360).  The  cavity 


motions  were  recorded  on  a  video  recorder  using  a  video  camera 


3.3  Test  Procedure 


The  range  of  test  parameters  considered  for  the  experiments  j 

were:  (a)  tunnel  pressure:  1  psia  -10  psia;  (b)  tunnel  velocity:  I 

i 

10  ft/sec  -45  ft/sec;  and  (c)  foil  oscillation  frequency:  0.2  Hz  [ 

-30  Hz.  Three  distinct  phases  of  investigations  were  undertaken  as  j 

described  below.  j 

I 

3.3.1  Video  Recording  of  Cavitv  Motion  ' 

The  unsteady  motions  of  the  cavity  were  recorded  on  a  video  tape 

i 

using  a  video  camera  for  several  freauencies  of  foil  oscillation  for  | 

I 

sufficiently  long  periods  to  record  many  complete  motion  cycles  at 

I 

each  frequency.  The  plan  view  of  the  cavity  was  recorded  and  it 
demonstrated  that  the  flow  was  nearlv  two  dimensional.  A  strobe  was 
used  to  obtain  tine  history  of  cavity  notion.  Care  was  taken  to  set 
the  frequency  of  the  strobe  at  a  value  that  was  different  from  the 
foil  oscillation  frequency  so  that  the  variations  of  cavity  length 
with  tine  could  be  determined.  In  this  way,  the  variation  of  cavity 
length  for  one  complete  cycle  of  oscillation  could  be  determined  by 
measuring  the  length  at  a  number  of  points  in  the  phase  of  the 
oscillation  as  explained  below, 

3.3.2  Source-Like  Nature  of  the  Oscillatlne  Cavitv 

These  tests  were  performed  by  measuring  the  unsteady  pressures 
at  three  field  points  at  the  tunnel  wall:  32.5  inches  upstream  and 
downstream  from  the  axis  of  oscillation  and  16.25  inches  upstream 
from  the  axi~  of  oscillation.  These  field  points  were  chosen  such 
that  any  "potential  flow"  pressures  at  these  wall  points  due  to  an 


i 


i 

I 


oscillating  cavity  being  present  in  the  cylindrical  test  section  were 
negligible  and  the  results  can  be  thought  of  as  far-field  free  stream 
fluctuations.  A  computer  code,  developed  by  Fernandez  (23),  was 
used  to  find  an  upper  bound  for  this  effect  assuming  the  shape  of  the 
cavity  to  be  an  ellipse.  Figure  3.6  shows  the  variation  of  pressure 
coefficient  along  the  length  of  the  tunnel  for  the  extreme  condition 
when  a  body  in  the  shape  of  an  ellipse  of  semi -major  axis  measuring 
3  inches  and  semi -minor  axis  measuring  1.5  inches  is  kept  in  the 
test  section.  At  each  field  point,  unsteady  pressures  were  measured 
for  several  tunnel  velocities  keeping  the  cavitation  number  constant. 

3.3.3  Measurement  of  Cavity  Pressures 

Steady  and  unsteady  cavity  pressures  were  measured  for  several 
foil  frequencies,  cavitation  numbers  and  static  angles  of  attack.  The 
steady  pressure  was  measured  from  the  pressure  tap  leading  into  the 
cavity.  This  pressure  tap  was  connected  to  the  pressure  transducer 
through  a  valve  venting  the  line  to  the  atmosphere  in  the  laboratory. 
By  closing  the  valve  only  after  the  formation  of  the  cavity,  the 
presence  of  water  bubbles  in  the  pressure  tap  line  was  completely 
eliminated.  Unsteady  pressure  data  were  analyzed  through  a  real  time 
analyzer. 

3.4  Test  Results  and  Discussion 

Sample  data  from  the  experiments  are  tabulated  in  Appendix  B.l. 
(A  complete  tabulation  of  the  experimental  data  is  available  upon 
reauest  from  the  author.)  The  cavitation  number  was  calculated 
based  on  the  measured  steady  cavity  pressure.  The  maximum  error 
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in  calculating  the  cavitation  number  is  estimated  to  be  4.8%  (see 
Appendix  B.3).  The  region  of  the  cavity  where  the  cavity  pressure 
was  measured  was  nearly  two-dimensional.  Hence  a  wall  effect 
correction  appropriate  to  such  flows  developed  by  Hsu  (24)  was  used 
to  correct  the  measured  cavitation  number.  An  example  of  this 
correction  procedure  is  contained  in  Appendix  B.4. 

3.4.1  Cavity  Length  History 

At  low  frequencies  of  foil  oscillation  (0.2  Hz  -  3.0  Hz)  it  was 
noticed  visually  that  the  cavity  behind  the  foil  seemed  to  follow  the 
motion  of  the  foil  without  any  noticeable  phase  difference  between 
input  foil  motion  and  the  resulting  cavity  length  caange.  However, 
as  the  driving  frequency  was  increased  beyond  5  Hz,  the  Dhase  lag  of 
the  cavity  motion  was  noticeable.  Beyond  about  20  Hz,  there  was  very 
little,  if  any,  motion  of  the  cavity  due  to  the  oscillation  of  the 
foil.  It  was  clear  from  these  visual  observations  that  the  cavity 
has  the  appearance  of  an  inertial  effect. 

The  cavity  length  history  was  obtained  by  playing  back  the  video 
tape  manually  on  a  large  television  screen  and  measuring  the  cavity 
lengths  at  equal  time  intervals,  the  time  interval  being  governed  by 
the  frequency  of  the  strobe.  The  magnification  of  the  cavity  on  the 
television  screen  was  properly  accounted  for  in  obtaining  true  cavity 
lengths.  The  cavity  length  measurement  on  the  screen  was  at  times 
hampered  by  the  "not  so  clear"  terminal  end  of  the  cavity.  However, 


all  the  measurements  were  consistently  made  and  the  cavity  lengths 
were  measured  within  an  accuracy  of  three  millimeters. 


The  tine  histories  of  cavity  motions  for  a  cavitation  number  of 
0.47  at  four  different  foil  frequencies  are  shown  in  Figures  3.7(a) 


through  3.7(d).  It  is  apparent  from  these  figures  that  the  motion 
of  the  cavity  is  far  from  being  sinusoidal.  The  Fourier  components 
of  these  time  histories  are  shown  in  Figures  3.8(a)  through  3.8(d), 
respectively.  It  is  observed  from  these  results  that  the  cavity 
motion  consists  of  frequencies  other  than  the  fundamental  driving 
frequency  and  in  particular  the  first  overtone  is  significant  in 
all  the  cases.  An  observation  of  the  spectra  at  different  foil 
oscillation  frequencies  shows  that  at  0.5  Hz,  the  motion  of  the 
cavity  is  close  to  harmonic  motion  whereas  at  1  Hz  the  great 
increase  in  the  first  overtone  of  the  unsteady  cavity  is  apparent. 

As  the  frequency  is  increased,  the  cavity  motion  is  again  close  to 
harmonic  motion  as  evidenced  from  the  spectra  in  Figures  3.8(c)  and 
3.8(d).  It  is  of  interest  to  compare  qualitatively  these  amplitudes 
of  harmonics  of  cavity  motion  with  the  linearized  unsteady  cavity 
theory.  If  we  denote  £(t)  as  the  cavity  length  at  an  instant,  we 
may  write  (see  Chapter  5) 

*(t)  -  ,  (3#1) 

where  a(t)  is  the  angle  of  attack  and  a  is  the  cavitation  number. 

For  the  oscillating  foil  motion,  we  may  write 

a(t)  =  oQ(l  +  ee'^)  ,  (3.2) 

where  ea^  denotes  the  amplitude  of  pitching  motion.  Recognizing 
the  fact  that  e  is  a  small  number,  we  can  write  for  £(t), 
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i(t)  *  (i  +  |  ejut  -  ej2uC  +  ...)  .  (3.3) 

The  similarity  between  the  measured  harmonic  contents  of  cavity 
length  and  those  obtained  from  Eq.  (3.3)  is  now  evident.  Figure  3.9 
shows  the  variation  of  amplitude  of  unsteady  cavity  motion  (defined 
as  the  peak  to  peak  amplitude  normalized  with  respect  to  steady 
cavity  length)  as  a  function  of  the  reduced  frequency.  It  is  clear 
that  the  amplitude  of  unsteady  cavity  motion  decreases  rapidly  with 
increasing  frequency.  It  is  also  of  interest  to  observe  that  the 
amplitude  of  unsteady  cavity  motion  peaks  at  a  particular  value  of 
reduced  frequency.  Several  attempts  were  made  to  measure  the  phase 
difference  between  the  cavity  motion  and  foil  motion.  3ecause  of  the 
extreme  difficulty  involved  in  such  measurements,  it  is  not  surprising 
that  no  meaningful  data  were  obtained,  using  the  present  measurement 
technique. 

3.4.2  Source-Like  Mature  of  Oscillating  Cavitv 

Figures  3.10  and  3.11  show  the  plots  of  the  relative  root  mean 
square  pressure  level  in  decibels  at  32.5  inches  upstream  of  the  pitch 
axis  as  a  function  of  the  test  section  velocity.  The  foil  oscillation 
frequency  is  relatively  low,  being  less  than  2  Hz.  It  is  seen  that 
the  unsteady  pressure  varies  as  the  sixth  power  of  the  test  section 
velocity.  This  indicates  that  the  oscillating  cavity  behaves  like  a 
dipole  source  even  at  these  low  frequencies  [see  for  example, 

Donald  Ross  (26)  and  Billet  and  Thompson  (27)].  Figure  3.12  shows  a 
similar  plot  of  unsteady  pressure  versus  test  section  velocity  at  a 
point  32.5  inches  downstream  of  the  pitch  axis.  Here  again,  the 
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dipole  behavior  of  the  oscillating  cavity  is  very  clear.  Figure  3.13 
depicts  the  measurements  at  a  location  16.325  inches  upstream  of  the 
pitch  axis.  At  an  oscillation  frequency  of  10  Hz,  the  unsteady 
pressure  varies  as  the  eighth  power  of  the  test  section  velocity. 

This  variation  suggests  that  the  oscillating  cavity  behaves  like  a 
quadrupole  at  this  frequency.  This  is  possibly  due  to  the  phase 
lag  that  exists  between  the  motion  of  the  foil  and  the  notion  of 
the  cavity  at  high  frequencies.  The  oscillating  foil  and  the 
oscillating  cavity  each  constitute  a  dipole  source  and  with  the 
adequate  phase  lag  that  exists  between  them  at  high  frequencies, 
they  jointly  constitute  a  quadrupole  source. 

3.4.3  Cavitv  Pressure  Measurements 

The  measurement  of  steady  cavity  pressures  as  well  as  the 
unsteady  cavity  pressures  were  accurate  to  within  27..  The  steady 
cavity  pressure  is  normalized  with  respect  to  the  ambient  vapor 
pressure  of  water.  The  variation  of  normalized  steady  cavity 
pressure  as  a  function  of  the  parameter  of  a  is  shown  in  Figure  3.14. 
The  average  value  of  the  dissolved  gas  content  in  the  tunnel  water 
was  9  ppm  by  moles  and  the  average  value  of  the  free  air  in  the 
tunnel  water  was  4  ppm  by  volume  during  these  tests.  The  steady 
cavity  pressure  is  seen  to  be  always  greater  than  the  vapor  pressure 
of  water  at  the  ambient  temperature.  It  is  observed  that  larger 
cavities  (lower  cavitation  numbers  and/or  higher  angles  of  attack) 
sustain  higher  cavity  pressures  than  smaller  cavities.  In  the 
range  of  of  a  between  2  and  6,  it  is  observed  that  the  steady  cavity 
pressure  has  a  nearly  linear  relation  with  a/a.  Figure  3.15  shows  the 
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variation  of  normalized  steady  cavity  pressure  at  a  steady  angle  of 
attack  of  8  degrees  as  a  function  of  the  reduced  frequency  parameter 
k'.  [Parkin  (28)  has  shown  that  for  cavitating  oscillating  foils, 
the  appropriate  non-dimensional  frequency  is  k'  =  k//l  +  <j,  where  k 
is  the  reduced  frequency  in  the  non-cavitating  case,  k  =  Ub/U  .]  The 

CO 

amplitude  of  oscil  ation  is  2  degrees  peak  to  peak.  It  is  observed 
that  the  steady  cavity  pressure  increases  with  frequency,  reaches  a 
maximum  value  and  then  decreases  with  further  increase  in  frequencv. 
This  can  possibly  be  explained  by  the  entrainment  process  behind  an 
oscillating  cavity.  As  the  oscillation  frequency  Increases,  it 
appears  as  though  the  re-entrant  jet  behind  the  cavitv  becomes 
increasingly  staoie  until  a  certain  irecuency  is  reac.oea.  The  staoie 
re-encranc  jec  reduces  air  entrainment  behind  the  cavity  ana  hence 
the  cavity  is  able  to  sustain  higher  cavity  pressures.  This 
situation  however  changes  at  high  frequencies  of  foil  oscillation 
wherein  the  re-entrant  jet  becomes  unstable  due  to  the  high  freauencv 
surface  waves  interacting  with  the  re-entrant  jet.  At  high  foil 
oscillation  frequencies,  the  cavity  remains  essentially  stationary  and 
therefore  the  entrainment  rate  is  equal  to  the  gas  diffusion  rate 
resulting  in  constant  cavity  pressure.  In  other  words,  so  far  as 
steady  cavity  pressure  is  concerned,  very  low  frequency  oscillations 
of  the  cavity  and  very  high  frequency  oscillations  of  the  cavity  have 
the  same  effect.  It  is  quite  interesting  to  observe  that  the 
frequency  at  which  the  maximum  steady  cavity  pressure  occurs  and  the 
freauency  at  which  the  maximum  change  in  cavity  length  occurs  (see 
Figure  3.9)  are  quite  nearly  the  same,  at  k'  = 
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The  variation  of  unsteady  cavity  pressure  normalized  with 
respect  to  the  steady  cavity  pressure  as  a  function  of  the  reduced 
frequency  is  shown  in  Figure  3.16.  Here  the  steady  angle  of  attack 
is  15  degrees  and  the  amplitude  of  foil  oscillation  is  4  degrees 
peak  to  peak.  It  is  clear  that  the  unsteady  pressure  peaks  at  a 
reduced  frequency  of  about  0.05  and  then  drops  off  rapidly  with 
increase  in  frequency.  This  behavior  of  unsteady  cavity  pressure 
with  frequency  may  be  explained  from  the  unsteady  cavity  length 
data.  At  very  high  frequencies,  there  is  very  little  motion  of  the 
cavity  and  hence  the  cavity  pressure  is  verv  small.  These  trends 
in  unsteady  cavity  pressures  acree  well  witn  r.cose  or  Jiang  ana 
Lee'ney  (17).  Figure  3.17  similarly  shows  the  normalized  unsteadv 
cavity  pressure  as  a  function  of  recucea  frequency  at  a  steady 
angle  of  attack  of  8  degrees  and  amplitude  of  oscillation  of 
2  degrees  peak  to  peak.  Here  the  peak  unsteady  cavity  pressure 
occurs  at  a  reduced  frequency  of  0.16.  This  is  to  be  expected 
because  the  cavity  at  a  foil  angle  of  3  degrees  is  smaller 
compared  to  the  cavity  at  15  degrees  and  hence  it  has  a  higher 
resonant  frequency. 

It  is  interesting  to  observe  from  Figure  3.14  that  the  linear 
relation  between  the  cavity  pressure  and  the  parameter  a/a  enables 
one  to  estimate  the  cavitation  number,  based  on  cavity  pressure,  in 
terms  of  free  stream  static  pressure,  the  vapor  pressure  at  the 
ambient  temperature,  the  free  stream  dynamic  pressure  and  the  angle  of 
attack: 
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3.4.4  Recapitulation 


These  experiments  on  unsteady  cavities  have  given  us  some  insight  into 
their  behavior.  The  variations  in  harmonic  content  for  the  motion  of 
the  cavity  have  been  clearly  established.  The  source-like  nature  of 
the  oscillating  cavity  has  been  characterized  as  a  dipole  source  at 
low  frequencies  of  oscillation  and  approaches  a  quadrupole  source  at 
high  frequencies.  One  would  expect  to  see  a  monopole  contribution  as 
well  because  of  the  cavity  volume  fluctuations  and  Fig.  3.13  suggests 
such  a  possibility.  However,  this  experimental  fragment  is 
overwhelmed  by  the  bulk  of  the  presently  available  data  favoring 
higher-order  singularities.  Since  the  cavitv  and  foil  general!'.’  are 
rather  large  for  a  six-inch  diameter  test  section,  the  influence  of 
tunnel  walls  on  these  experimental  findings  remains  to  be 
investigated.  Low  cavitation  numbers  and/or  high  angles  of  attack 
result  in  higher  cavity  pressures.  Unsteady  cavity  pressure  is  seen 
to  be  significant  only  at  low  frequencies  of  oscillation. 


CHAPTER  4 


TWO-DIMENSIONAL  UNSTEADY  DIFFUSION  ACROSS 
A  PERMEABLE  CAVITY  SURFACE 

4.1  Introduction 

From  the  previous  studies  on  cavity  pressure  fluctuations  by  Jiang 
and  Leehey  (17)  and  the  present  experimental  investigations  described  in 
Chapter  3,  it  is  clear  that  gaseous  diffusion  plays  a  very  important 
role  in  the  behavior  of  unsteady  cavities.  Parkin  (8)  studied  steady 
diffusion  across  a  two-dimensional  cavity  surface,  while  Brennen  (9) 
considered  steady  diffusion  across  a  cavity  surface  from  a  turbulent 
flow. 

It  aooears  however,  that  a  corresooncing  analysis  for  unsteacv 
diffusion  in  an  unsteadv  cavity  flow  is,  to  the  best  of  our  knowledge, 
non-existent.  Unsteady  gaseous  diffusion  into  a  cavity  may  result  from 
two  distinct  mechanisms.  First,  since  the  cavity  length  fluctuates, 
the  gaseous  diffusion  into  the  cavity  fluctuates.  Second,  since  the 
cavity  gas  pressure  itself  fluctuates,  the  concentration  gradient 
across  the  cavity  surface  is  a  function  of  time.  In  this  chapter, 
we  formulate  and  solve  analytically,  the  unsteady  gaseous  diffusion 
across  a  two-dimensional  unsteady  cavity  surface.  The  results  of 
this  analysis  will  give  a  better  understanding  of  the  parameters 
that  affect  the  unsteady  gaseous  diffusion  across  a  cavity  surface. 

4 . 2  Assumptions 

It  is  assumed  in  the  following  analysis  that  the  flow  is  in¬ 
compressible,  two-dimensional,  isotropic  and  turbulent.  (Turbulence 
is  inherently  three-dimensional  in  nature.  Nonetheless,  a  two- 
dimensional  flow  field  with  unit  depth  may  be  considered  for  diffusion 
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analysis.)  Boussinesq's  hypothesis  (29)  is  used  in  turbulence  modeling 
and  G.  I.  Taylor's  (30)  statistical  theory  of  turbulence  is  used  in 
estimating  the  mass  diffusivity  for  the  flow  field.  It  is  assumed  that 
the  mechanism  for  turbulent  momentum  transfer  and  turbulent  concentration 
transfer  are  identical.  It  is  also  assumed  that  the  gaseous  diffusion 
occurs  instantaneously  across  the  cavity  surface. 


4.3  Theoretical  Analysis 

Gaseous  diffusion  will  occur  when  there  exists  a  dissolved -gas 
concentration  difference  between  the  free  stream  and  the  liquid  on  the 
cavity  surface.  If  at£  denotes  the  dissolved  gas  content  in  the  free 
stream  in  parts  per  million,  by  moles,  then  by  Henrv 's  law,  the  maximum 
partial  pressure  of  gas  in  the  free  stream  P_^  is 

PFS  -  V  •  <1-1 

where  3  is  the  Henry's  law  constant.  If  p_(t)  is  the  instantaneous 

b 

partial  pressure  of  non-condensable  gas  in  the  cavity,  then  the  mean 
concentration  difference  c(t)  expressed  in  moles  is 


:<t)  -  (a,  -f)  . 


We  assume,  as  is  customary,  that  the  mechanism  for  turbulent 
diffusion  of  gas  in  the  liquid  and  across  the  wall  and  the  turbulent 
momentum  transfer  are  similar  and  that  gradients  in  mean  velocity  can 
be  neglected.  Then,  the  turbulent  diffusion  is  due  solely  to  the 
gradient  of  the  mean  concentration  c(x,y).  Consider  a  turbulent  flow 
field  having  a  uniform  mean  velocity  in  the  positive  x-direction,  as 
shown  in  Figure  4.1.  Let  the  instantaneous  concentration  at  a  point 
(x,y)  be  C(x,y,t).  Also,  let  C  be  the  concentration  at  large  distances 
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from  the  origin  and  be  the  saturation  concentration  on  the  cavity 

wall.  Let  c  *  C  -  C  be  the  concentration  difference.  From  Fick's  law 

00 

of  diffusion  and  conservation  of  mass  [see,  for  example,  Batchelor  (31)] 
it  follows  that 

(U  •  V)c  +  =  kV^c  (4.3) 

2 

where  V  denotes  the  Cartesian  Laplacian  operator  and  k  is  the  molecular 
diffusivity. 

Let  the  instantaneous  velocity  components  be  represented  by 
(U  +  u'),  v'  and  w'  in  the  three  orthogonal  directions  and  let  c  =  c  +  c 
where  c  is  the  time  mean  concentration  difference.  The  comDonents  u',  v ' 
and  w'  are  the  instantaneous  velocity  fluctuations  and  c'  is  the 
instantaneous  concentration  difference.  Then,  Eq .  (4.3)  becomes  [see, 
for  example,  Goldstein  (32)] 
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+  u  - - <7  c  +  — — 

3t  ®  ox  3t 


U 

a 


3c' 

3x 


+  u 


3c 

3x 


+  u 


3c' 

3x 


+  v  ' 


3£ 

3y 


aC 

3y 
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aZ  dz  J 


(4.4) 


Time-averaging  the  above  equation,  we  obtain 


3c 


3c  tT  3c  3  (  3c  — ; — 3 r  3c  — ; — 3 r  3c  — ; — p, 

3t  ®  3x  3x  ^  3x  '  3y  v  3y  '  ^  j 


3z  ^  3z 


(4.5) 

where  the  overbar  indicates  time  averaging.  Equation  (4.5)  shows  that  in 


turbulent  flows  the  eddy  concentration  transport  terms  -  u'c',  -  v'c'  and 


3c  3c 


-  w'c'  add  to  the  molecular  concentration  transports  ic  — ,  k 


and 


<  —  respectively.  The  eddy  transport  terms  are  large  compared  to 
3z 

molecular  transport  terms  (32).  Therefore,  the  latter  are  neglected  and 
Eq.  (4.5)  reduces  to 


i£.+  u  !£■ 
at  ®  ax 


■ir(u'c')  v'c')  ~!r(-w,c') 


(4.6) 


We  now  restrict  our  study  to  two-dimensional  flows.  We  represent 
the  turbulent  shear  stress  pu'c'  and  pv'c'  using  Boussinesq's  hypothesis 


r  ,  ou 

-  pu  c  =  ^  —  , 


(4.7a) 


-  pv'c'  =  A,  —  , 

J.  oV 


(4.7b) 


where  A^,  is  the  exchange  coefficient  or  eddy  mass  conductivity.  Let  vt 
denote  the  eddy  mass  diffusivity.  Then, 


vt  =  p  , 


(4.8) 


where  p  is  the  density  of  the  liquid.  Substitution  of  Eqs.  (4.7)  and 
(4.8)  in  (4.6)  yields 


—  —  2—2— 

3c  8c  ,  a  c  a  c . 
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(4.9) 


The  boundary  conditions  pertinent  to  the  two-dimensional  unsteady 
diffusion  are: 


=  0  ;  y  =  0,  x  <  0  ,  x  >  i 


(4.10a) 
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c  =  c  (1  +  c(t))  ;y  =  0,0<x<£ 
o 


(4.10b) 


where  £(t)  is  the  unsteady  cavity  length  given  by 


£(t)  =  £  (1  +  e£(t ) )  , 

o 


(4.11) 


where  e  is  a  small  quantity. 


We  now  let 


c(x,y,t)  =  cg(x,y)  +  cu(x,y,t)  , 


(4.12) 


where  cs  and  cu  represent  respectively,  the  steady  part  and  unsteady  part 


of  the  concentration  c. 

Substitution  of  Ea .  (4.12)  in  Ea .  (4.9)  results  in  two  eauations 
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(4.13b) 


The  steady  partial  differential  equation  (4.13a)  satisfies  the  steady 
part  of  the  boundary  condition  (4.10b)  viz., 


c  =c  ;  y  =  0,  0<x<£ 

so  o 


(4.14) 


The  unsteady  partial  differential  equation  (4.13b)  satisfies  the  unsteady 
part  of  the  boundary  condition  (4.10b)  viz., 


v.v! 
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=  CQc(t)  ;  y  =  0,  0  <  x  <  £(t) 


(4.15) 


The  steady  partial  differential  equation  has  been  solved  by  Parkin  (8). 

The  solution  of  Eq.  (4.13b)  subject  to  the  boundary  conditions  (4.15)  is 
obtained  in  stages.  First,  a  point  source  solution  that  has  the  proper 
radial  symmetry  and  behavior  is  obtained.  Second,  the  point  source 
solution  is  used  to  formulate  an  integral  equation  relating  the 
concentration  at  a  fixed  point  and  the  mass  flux  per  unit  length  along  the 
cavity.  Third,  the  concentration  difference  at  the  cavity  surface  is  used 
to  solve  the  integral  equation. 

The  fundamental  solution  of  Eq .  (4.13b)  for  an  oscillating  unit 
source  at  the  origin  is  (See  Appendix  C.l) 


c(3t-y:!0  -tis; 


K  (R  Vx“  -  v"  )  e 
o 


(E-jk)x 


(-.16) 


where  K^z)  is  the  modified  Bessel's  function  of  the  second  kind.  It 
should  be  noted  that  r,  x  and  y  appearing  in  Equation  (4.16)  are 


dimensionless  and  that  R  =  - —  ,  a  "diffusive  Reynold's  number." 

2vt 

We  now  introduce  the  harmonic  oscillations  of  the  unsteady  boundary 
conditions  resulting  from  the  angle  of  attack  oscillations  at  a  circular 


frequency  of  m  radians  per  second.  Suppose  that  the  cavity  length 
oscillations  lag  the  foil  oscillations  (see  Chapter  5)  by  a  phase  angle  of 
4).  Then,  the  unsteady  cavity  length  at  x  =  £0  is  given  by  e£0  ei^ajC+'t>)  or 
t  lQ  ei  $  el . 


m 
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The  gas  pressure  within  the  cavity  Pq(0  may  be  written  as 


PG(C)  =  PGo(1  +  5  P<'t)')  ’ 


(4.17) 


where  <5  is  a  small  quantity. 

With  the  assumption  that  the  harmonic  cavity  pressure  fluctuations 
lag  the  foil  oscillations  by  $a,  we  may  write  from  Eq.  (4.2), 


c(t)  =  (a  -  a  )  -  a  <5p  , 

1  a  a 


(4.18) 


So  that  the  concentration  fluctuation  cu(t)  across  the  permeable  surface 


is  given  by 


i  .  e .  , 


c  (t )  =  - 


’  ( it  +  0  ) 


, .  V  -  -  (  -T  -  2  )  jut 

:(:)=-  a  c  e“  a  e- 

u  a 


(4.19) 


where  <ia  is  the  phase  angle  between  foil  motion  and  cavity  pressure 
fluctuations.  Equation  (4.19)  clearly  indicates  a  phase  lag  of  tt  radians 
between  the  concentration  fluctuations  and  cavity  pressure  fluctuations. 
The  steady  concentration  difference  cQ  at  the  permeable  surface  is 


c  =  (a,  -  a  ) 
o  la 


(4.20) 


Let  f(x;k)  be  the  mass  flux  per  unit  length.  Then  the  function  c(x,y;k) 
can  be  expressed  in  terms  of  f(x;k)  as  follows: 


;(x,y;k)  =  y-— -  /  f  ( C;k)e  (R  jk)  (x'°K  (R  (x-O2  +  y 2  )  d£  , 
2ttAT  o  0 


(4.21) 


where 


--2-(l  +  ceJ*)  . 


If  we  now  require  that  c(x,o;k)  =  cQ  -aaS  eJ  ^  in  the  interval  0  <  x  <  £., 
the  source  strength  is  determined  from  the  integral  equation 


2  m 
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+  /  fU;k)e"R(5  x)  KQ(R(5-x))ejk(?'x)d5]  ,  (4.22) 

x 

where  the  positive  branch  of  the  square  root  has  been  taken.  The  function 
c(x,y;k)  also  satisfies  the  conditions  that  c+o  as  (x2  +  y2)  +  <»  and  that 

3c/3y(x,o)  =  0  when  x  <  0  and  x  >  l  . 

In  Eq .  (-.22)  we  have  combined  the  steaay  as  well  as  the  unsteady 
boundary  conditions  for  both  cavity  length  and  cavity  pressure.  The 
reason  for  this  is  that  the  point  source  solution  in  the  unsteady  case 
reduces  to  the  steady  source  solution  as  the  reduced  freauencv  k  goes  to 
zero.  This  fact  also  implies  that  5  =  5(k,o)  ,  e  =  e(k,o)  and  that  5  and 
s  both  vanish  as  k  *  o. 

In  order  to  determine  f(x;k),  we  shall  make  use  of  the  fact  that  R  is 
a  large  number  for  the  present  study.  Therefore,  Kq  in  Eq .  (4.22)  can  be 
replaced  by  the  first  term  of  Its  asymptotic  expansion, 

K0(z)  =  e~Vi7  '  (4.23) 

A  further  simplification  can  be  obtained  in  the  second  integral  in  Eq . 
(4.22)  by  noting  that  the  strong  negative  exponential  will  cause  f(£;k)  to 
contribute  to  the  integration  only  near  £  =>  x.  Therefore,  we  replace  the 

second  integral  in  Eq .  (4.22)  by  the  approximate  value  f  (x;k)  — —  . 

2  /  2R-jk 
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With  these  approximations,  Eq .  (4.22)  may  be  written  as 
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If  the  Laplace  transform  of  f(x;k)  is  denoted  by  F(s;k),  we  can  transform 


Eq.  (4.24)  to  find 


F(s;k)  = 
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(4.25a) 


Rearranging  terms,  we  have 
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The  inverse  of  this  transform  gives  (see  Appendix  C.2)  for  f(x;k')  the 


result , 


2  —  — 

f(x;k)  =  0[S-  -  e  ^  erfc(b/x)  -  erf/ax  ]  ,  (4.26) 
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(4.27a) 


a  =  jk  , 


(4.27b) 


b  =  /2R  -  jk  , 


(4.27c) 


a  =  2jk  -  2R  , 


(4. 27d ) 


l  =  £(t)/b  as  before, 
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and  erf(x)  and  erfc(x)  denote  the  error  function  and  complimentary  error 


function  respectively.  We  note  here  the  fact  that  the  function  f(x;k) 


which  represents  the  mass  flux  per  unit  length  is  independent  of  the 


cavity  length. 


The  rate  at  which  the  mass  is  diffused  per  unit  width  along  the 


entire  length  of  the  permeable  interval  is  obtained  from 


/  f (x,k)  dx  , 


(4.28) 


.  I  ,2  __  _ 

i.e.,  M(k)  -  Q  f  (7  -  —  e  201  erfc  b/x —  erf  /ax  )  dx  . 


(4.29) 
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The  integral  in  Eq .  (4.29)  has  been  evaluated  in  Appendix  C.3  and  the 


result  is 
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(4.30) 


Equation  (4.30)  gives  the  instantaneous  unsteady  mass  flow  rate  along  the 


cavity  surface.  We  now  compare  this  instantaneous  mass  flow  rate  relative 


to  Mq,  the  mass  flow  rate  obtained  in  the  steady  limit.  By  letting 


the  reduced  frequency  k  go  to  zero,  we  obtain  the  steady  limit  for 


mass  flow  rate  as 
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We  can  now  write  the  instantaneous  mass  flow  rate  in  terms  of  the  quasi¬ 
steady  mass  flow  rate  as 


M  *  M 
o 


[1  +  r(k)e 


(4.32) 


where  Mgr(k)  represents  the  absolute  value  of  the  frequency  dependent  mass 
flow  rate  and  ij>  is  the  phase  angle  by  which  the  instantaneous  mass  flow  rate 
leads  or  lags  the  steady  value. 

It  is  of  interest  now  to  look  at  the  relative  effects  of  varying  the 
concentration  and  cavity  length  on  the  unsteady  mass  flow  rate.  3y  a 
simple  sensitivity  analysis  of  the  auasi  steady  mass  flow  rate  (see 
Appendix  D.l)  we  obtain 


AM  ^  ^o  +  A* 

M  Co  2l 
o  o 


(4.33) 


where  the  A  before  a  quantity  represents  a  small  change  in  that  particular 

quantity.  For  the  experimental  investigation  reported  in  Chapter  3,  the 
Ac 

maximum  value  of  -  —  is  of  the  order  of  5%,  whereas  the  maximum  value  of 
c 

o 

A 

■7—  is  of  the  order  of  20 7..  Thus,  it  is  clear  that  for  the  experiments 
4. 

o 

reported  in  Chapter  3,  the  contribution  to  gas  diffusion  across  the  cavity 
from  changing  cavity  lengths  and  changing  cavity  gas  pressure  are  equally 
significant. 


4.4  Recapitulation 


We  have  developed  an  analytical  model  for  the  prediction  of  gaseous 
diffusion  across  a  two-dimensional  unsteady  cavity  surface.  This  model 
takes  into  account  the  change  in  cavity  length  as  well  as  changes  in 
cavity  pressure  in  predicting  the  gas  diffusion  across  the  cavity  surface. 
The  expressions  for  mass  flux  per  unit  length  and  mass  diffusion  rate 
across  the  entire  cavity  length  reduce  to  those  obtained  by  Parkin  (8) 
when  the  reduced  frequency  k  is  zero.  The  instantaneous  mass  flow  rate 
across  the  unsteady  cavity  may  be  easily  obtained  from  the  auasi-steady 
mass  flow  rate  and  the  lag  function.  We  reiterate  here  the  general  fact 
that  this  unsteady  gaseous  diffusion  analysis  is  valid  onlv  for  harmonic 
variations  in  cavity  gas  pressure  and  possibly  in  cavity  length. 
Nonetheless,  this  analysis  could  be  extended  to  encompass  a  general 
nonharmonic  motion  of  the  cavity  by  Fourier  representation  of  the  cavity 
length  and  cavity  gas  pressure  fluctuations,  provided  that  one  first 
accounts  for  the  nonsteady  cavity  closure  conditions  to  be  considered 


CHAPTER  5 


A  STUDY  OF  NONSTEADY  CAVITY  CLOSURE 

5.1  Introduction 

The  analysis  of  a  two-dimensional  unsteady  cavity  of  finite 
extent  (i.e.,  non-zero  cavitation  number)  has  been  studied  among 
others,  by  Parkin  (28),  Wu  (5),  Guerst  (33)  and  Acosta  and 
Furuya  (34).  The  fact  that  a  cavity  volume  change  in  two- 
dimensional  unbounded  flow  would  create  infinite  pressures 
puzzled  many  researchers  and  was  thought  to  be  non-physical  at 
the  time.  Wu  (5)  points  out  the  fact  that  such  cavity  volume 
changes  in  two-dimensional  flow  are  indeed  physically  plausible 
and  mathematically  they  represent  ''inner"  flows  of  some  three- 
dimensional  outer  flows.  Acosta  and  Furuya  (34)  rightly  exemplify 
this  fact  by  the  existence  of  unsteady  cavities  in  water  tunnel 
experiments. 

The  presence  of  gases  In  an  unsteady  cavity  flow  makes  the 
study  of  such  flows  even  more  interesting  and  introduces  more 
complexities.  Recently,  Kato  (35)  has  pointed  out  the  general  lack 
of  knowledge  that  exists  in  understanding  transport  phenomena  at 
the  cavity  surface.  The  results  of  the  unsteady  cavity  flow 
experiments  reported  in  Chapter  3  have  clearly  shown  that  the  cavity 
pressure  does  not  remain  constant  during  the  oscillation  of  a  finite 
cavity  and  indeed  even  the  steady  gas  pressure  of  an  oscillating 
cavity  depends  on  the  frequency  of  oscillation!  Clearly,  the 
unsteady  gas  diffusion  plays  a  role  in  the  dynamics  of  an  unsteady 
cavity. 


The  idealized  closure  condicion  in  cavity  flow  means  that  at 
every  instant,  the  cavity-body  system  is  a  closed  body.  Several 
models  have  been  postulated  in  the  past  for  the  closure  of  an 
unsteady  cavity.  Among  then  are  the  kinematically  closed  model 
proposed  by  Tulin  (36),  the  constant  volume  model  proposed  by  Guerst 
(33),  the  fixed  cavity  terminus  model  proposed  by  Leehey  (37)  and 
the  generalized  cavity  closure  model  proposed  by  Parkin  (28). 

Acosta  and  Furuya  (34)  have  demonstrated  that  Parkin's  model  which 
involves  the  evaluation  of  a  difficult  retarded  integral  has  the 
proper  quasi-steady  limit  when  the  reduced  frequency  goes  to  zero. 
The  closure  condition  postulated  by  Parkin  has  been  and  remains  to 
be  difficult  to  evaluate  from  first  principles  of  the  analysis 
although  Parkin  proposed  a  step-by-step  numerical  scheme  which  has 
never  been  employed.  In  the  present  study,  we  shall  evaluate  the 
phase  angle  between  foil  motion  and  cavity  motion  by  means  of  a 
simplified  kinematic  nonsteady  closure  condition  and  then  relate  the 
cavity  pressure  fluctuations  and  cavity  volume  fluctuations.  Such  an 
analysis  should  lead  to  (i)  a  better  understanding  of  the  gaseous 
entrainment  behind  nonsteady  cavities  and  (ii)  relate  observed 
pressure  pulsations  to  unsteady  gaseous  entrainment. 

5.2  Steady  Cavity  Closure 

We  first  look  at  the  steady  cavity  closure  when  there  exists  a 
certain  amount  of  gas  in  addition  to  water  vapor  within  the  cavity. 
Continuity  of  mass  flux  of  gas  through  the  cavity  dictates  in 
general  that  the  mass  flow  rate  of  gases  through  the  cavity  must 
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equal  Che  rate  of  entrainment  behind  the  cavity  plus  a  possible  rate 
of  change  of  air  mass  inside  the  cavity.  In  the  limit  of  steady 
flow  this  last  part  will  vanish.  Thus,  we  can  estimate  the  rate  of 
entrainment  of  gases  behind  the  cavity  by  estimating  the  gas 
diffusion  rate  through  the  cavity  surface. 

For  a  flat  plate  at  a  small  angle  of  attack  ct^,  we  can 
estimate  the  cavity  length  at  a  cavitation  number  c  (28)  as 
(see  Appendix  E.l) 


(5.1) 


From  the  steady  cavity  pressure  measurements  reported  in  Chaoter  3 
(see  Figure  3.14),  we  can  approximate  the  dependence  or  cavity  gas 
pressure  and  a/ as 

Pr  - 

-=■  •-  0.0976  —  +  0.74  .  (5.2) 

P  ct- 

rv  0 

If  we  denote  by  a^,  the  amount  of  dissolved  gas  in  tunnel  water, 
the  concentration  difference  (in  ppm)  that  exists  between  the 
free  stream  and  cavity  surface  may  be  expressed  as 

e  “  Oj  -  Pq/B  .  (2.9) 


Thus  we  observe  from  Eqs.  (5.2)  and  (2.9)  that  the  concentration 
difference  c  is  a  function  of  a/a^.  We  can  obtain  the  mass  flux 
into  the  cavity  from  the  upper  surface  as  well  as  the  lower  surface 
of  the  two-dimensional  cavity  from  a  modification  of  Parkin's 
diffusion  model  as 
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f 


(5.3) 


where  vt  is  the  turbulent  eddy  mass  diffusivity  and  where  £g 
represents  the  cavity  length  on  the  upper  surface  and  represents 
the  cavity  length  on  the  lower  surface.  Figure  5.1  shows  the 
comparison  between  the  gas  entrainment  rates  when  the  cavity  gas 
pressure  is  nonzero  and  for  a  hypothetical  purely  vaporous  cavity. 
This  fictitious  vaporous  cavity  is  used  only  to  bound. the  diffusion 
rate.  We  observe  that  the  entrainment  rate  decreases  with  the 
presence  of  gases  within  the  cavity  as  expected.  For  small  values 
of  or/ag,  the  difference  in  entrainment  rate  is  hardly  noticeable 
because  at  low  values  of  o/ag  the  mass  diffusion  rate  due  to  an 
increase  in  cavity  length  is  more  significant  than  the  decrease  in 
mass  flow  rate  resulting  from  increased  cavity  gas  pressure. 

5.3  Nonsteadv  Cavity  Closure,  Phase  Relations 

The  difficulties  encountered  in  measuring  experimentally  the 
phase  angle  between  the  foil  motion  and  cavity  motion  has  been 
alluded  to  in  Chapter  3.  Because  of  the  lack,  of  knowledge  that 
exists  in  the  literature  in  the  determination  of  the  present  phase 
relations,  we  turn  to  physically  motivated  arguments  aimed  at 
finding  the  phase  difference  between  cavity  length  oscillations  and 
hydrofoil  pitching  oscillations. 

Consider  a  two-dimensional  cavity  behind  a  foil,  as  shown  in 
Figure  5.2a.  The  foil  has  a  steady  angle  of  attack  aQ.  We  suppose 
that  the  foil  oscillates  in  pitch  about  B  with  a  frequency  of  u 
radians  per  second,  the  amplitude  of  oscillation  being  A«.  The 
angle  of  attack  «(t)  at  any  instant  is  then  given  by  «(t)  -  «o  +  A« 
sin  ut.  We  now  approximate  the  upper  surface  and  lower  surface  of 


the  cavity  by  the  slits  AC  and  BC  respectively  as  shown  in  Figure  5.2b. 
The  closure  condition  of  the  nonsteady  cavity  depends  on  kinematic 
waves  (28),  having  constant  celerity  qc  but  different  travel  distances 
on  the  upper  and  lower  surfaces,  terminating  at  a  single  point  C  behind 
the  cavity.  For  kinematic  waves  originating  simultaneously  from  points 
A  and  B  to  reach  C  at  the  same  instant,  there  must  exist  a  phase  lag 
$  in  the  wave  originating  from  the  lower  surface  relative  to  the  wave 
originating  from  the  upper  surface.  Since  a(t)  is  responsible  for 
these  waves,  we  can  say  that  the  phase  of  the  cavity  terminus 
oscillation  lags  that  of  the  driver  a(t). 

If  the  travel  distance  on  the  upper  surface  is  and  the  travel 
time  from  A  to  C  is  tu,  then  tu  is  given  by  iy/qj..  The  travel  time  on 
the  lower  surface  of  distance,  Z^,  is  given  by  t^  =  (Zy  -  b)/qc  if  the 
cavity  is  closed.  But  the  difference  in  travel  time  between  the  upper 
and  che  lower  surfaces  is  Cu  -  C«  *  At  *  ~  and  then  oj  At  *  —  . 

*  q_  q„ 

c  c 

Therefore  the  closed-cavity  termination  phase  lag  is  determined  by 

u  At  =  ■  —  .  (5.4) 

qc 

We  now  recognize  the  quantity  —  as  the  modified  reduced  frequency  k' 

qc 

and  write 

+  *  k'  .  (5.5) 

Thus  we  observe,  within  the  limitations  of  the  present  assumptions  that 
if  the  kinematic  waves  (28)  from  the  upper  surface  and  lower  surface 
terminate  at  one  point  at  the  same  instant,  the  phase  angle  equals  the 
modified  reduced  frequency.  This  result  conforms  with  the 
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qualitative  visual  observations  made  during  the  experimental 
investigation  viz.,  at  low  frequencies  there  is  hardly  any  phase 
difference  discernible  whereas  at  higher  frequncies,  there  seems  to  be 
a  noticeable  phase  difference  between  the  foil  motion  and  cavity 
motion. 

We  note  in  passing,  that  this  simple  analysis  is  consistent  with 
the  linearized  approach  of  Parkin  (28),  but  we  have  used  the 
preceeding  physical  argument  in  order  to  avoid  the  nasty  retarded 
closure  integral  introduced  in  (28).  If  we  were  to  relax  this 
consistency,  we  would  replace  the  factor  b  by  b  cos  «0,  from  which  we 
would  have  the  phase  angle  as  k'  cos  *0  instead  of  k'.  As  to  be 
expected,  in  the  limit,  for  small  oscillations  of  a  vertical  plate 
there  would  be  little  or  no  phase  difference  between  the  upper  and 
lower  surfaces  of  the  cavity. 

An  additional  aspect  of  the  present  approximate  argument  is  that 
the  value  of  At  is  independent  of  the  fundamental  frequency  of  the 
foil  oscillation,  the  important  factor  being  the  constant  value  of  qc. 
Because  of  this  fact  this  same  phase  lag,  $,  between  profile  and 
cavity  terminus  oscillations  oust  apply  to  the  higher  harmonics  of  the 
cavity  length  fluctuations  as  reported  in  Figure  3.8  above  and  as  also 
suggested  by  Eq.  (3.3).  this  means  that  this  phase  lag  angle  depends 
on  the  basic  foil  oscillation.  Cavity  oscillation  overtones  of  2w  and 
3u,  etc.,  would  not  have  phase  lags  of  2$,  3<J>,  etc,  as  might  be 
supposed  at  first  glance. 

5«A  Nonsteady  Cavity  Closure,  Entrainment  Rates 

If  the  cavity  ordinate  at  any  instant  is  r((x,t),  then  the  cavity 
closure  condition  is  expressed  as  (28), 


0 

:«• 


/  dii  *  0  Ct  fixed) 
body 


(5.6) 


Now,  we  shall  assume  that  the  cavity  behind  a  flat  plate  of  chord 
length  b,  at  an  angle  of  attack  a(t),  at  any  instant  t,  may  be 
approximated  by  a  part  of  an  ellipse  fitted  to  the  cavity  termination 
point,  £(t),  and  to  the  leading  edge  and  trailing  edge  separation 
points,  see  Figure  5.2a.  We  also  assume  that  the  length  of  the  cavity 
£( t)  is  measured  from  the  nose  of  the  flat  plate  as  shown.  Then,  we 
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take  the  equation  of  the  ellipse  to  be 


(x  -  h)2  (y  ~  k)2 

2  .2  * 

a  b 


(5.7) 


where  (h,k)  represents  the  center  of  the  ellipse  and  the  quantities  a 
and  b  represent  the  semi -major  axis  and  the  semi -minor  axis  of  the 
ellipse  respectively.  At  any  instant,  the  ellipse  will  pass  through 
the  points  (0,0),  (£,  -  y  sin  a)  and  (b  cos  a,  -  b  sin  a).  These 
three  conditions  lead  to 
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For  small  values  of  the  angle  of  attack  a,  we  can  write  from  Equations 


5.8a,  b  and  c. 


k  -  ~  f  b  , 


.  b 
h  *  j  , 


a  -  y  (1  +  4i(2.-l) ) 1/2 


.  _  ab  ,1  +  4£(£-l)  ' 

°  4  '  -  -  ' 

l(l  -  1) 


(5.9) 


(5.10) 


(5.11) 


(5.12) 


It  should  be  noted  that  the  cavity  length  appearing  in  Eqs.  (5.11)  and 
(5.12)  is  dimensionless,  (i/b).  From  Eqs.  (5.9)  and  (5.10)  it  is 
clear  that  for  small  angles  of  attack,  the  center  of  the  elliptical 
cavity  is  fixed  to  the  center  of  the  foil.  From  Fig.  5.2a  it  is 
apparent  that  the  cavity  volume  V  per  unit  depth  is  then  simply  the 
the  area  of  the  semi  ellipse 


where 


v  -  l  (1  -  -4)  , 

21 


a  =  aQ  (1  +  da  e^ut)  , 


t  -  £  (1  +  e  e^C  ej*)  . 
o 


(5.13a) 


Equation  (5.13a)  can  be  written  as 


ira  b2  _  .  (ee^  +  Aa  -  Aa/2£  )e^ut 

V- -  (l  --)[!+ - = - —  ].  (5.13b) 
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The  steady  cavity  volume  V0  from  Eq.  (5.13b)  is 


—2 
ira  b“ 
o 


P.-T)  • 


(5.13c) 


so  that  the  cavity  volume  V  at  any  instant  can  be  written  as 


V  -  Vq(1  +  veJai  )  , 


where 


(5.14) 


»  -  ♦  A«  -  42-)/(l  -4-)  . 


(5.15) 


We  now  look  at  the  conservation  of  mass  flux  of  gas  through  the 
control  volume  consisting  of  the  cavity  surface,  the  gas  within  the 
cavity  and  the  mass  entrained  from  the  cavity  terminus.  Let  m  be 
the  total  mass  of  gas  within  the  cavity  at  any  instant,  be  the 
mass  flux  of  gas  diffusing  into  the  cavity  through  its  surface  and 
E  be  the  mass  flux  of  gas  entrained  out  of  the  cavity  through  the 
cavity  end.  Then  we  can  write  the  continuity  of  mass  as 
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(5.16a) 


Let  represent  the  steady  gaseous  diffusion  rate  across  the 

o 

entire  cavity  surface.  Then,  from  Eq.  (5.16a),  the  gaseous  mass 


entrainment  rate  E  relative  to  the  steady  gaseous  mass  entrainment 


rate  is  given  by 
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(5.16b) 


By  using  a  quasisteady  approach,  the  incremental  mass  diffusion 
rate  due  to  cavity  pressure  changes  and  cavity  length  changes  may  be 
written  from  Eq.  (5.3)  as, 


(5.17) 


12 1  la 
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where  l  and  1  represent  the  steady  dimensionless  cavity 
uo  o 

lengths  on  the  upper  surface  and  lower  surface  respectively.  The 


quantity 


has  been  evaluated  in  Appendix  E.2  in  terms  of  a/a 


and  Aa  /a  and  the  result  is 
o  o 


-0.0976  (a/a  )(Aa  /a  )(p  /8) 

_ o _ o  O  V 

Pv  a 

a.-  (0.74  -  0.0976  —  ) 
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(5.18) 


Assuming  that  the  gas  within  the  cavity  behaves  isothermally,  one  may 
use  the  equation  of  state  (see  Appendix  E.3)  to  obtain  for  3m/3t, 
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where  R  -  specific  gas  constant  (-1718  ft. lb. /slug  #R), 

4  -  phase  angle  between  cavity  pressure  fluctuation  and 

foil  oscillation, 

T  -  absolute  temperature  of  the  gas  within  the  cavity. 

From  Eq.  (5.3)  Che  steady  gaseous  diffusion  rate  is  given  by 

Md  -  — —  [  V2R  iu  +  V2R  l  ]  ,  (5.20) 

O  V  IT  O  O 

where  R  in  Eq.  (5.20)  is  the  diffusive  Reynolds'  number, 

R  -  U  b/2v  .  Eqs.  (5.17)  through  (5.20)  may  now  be  substituted  in 

•  t 

Eq.  (5.16b)  to  obtain  the  total  entrainment  rate  E/Md  as 
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A  polar  plot  of  the  total  entrainment  rate  for  four  different 
reduced  frequencies  is  shown  in  Figure  5.3.  The  experimental  data 
points  from  Figure  3.8  and  3.14  on  cavity  length  amplitudes  and  a 
quasisteady  gas  pressure  amplitude  based  on  the  linear  regression  of 
Figure  3.14  have  been  used  in  obtaining  these  calculated  data  points. 
(See  Appendix  E.3  and  E.4  for  sample  calculations  and  tabulated 
results).  It  oust  be  noted  that  the  unsteady  entrainment  rate 
amplitudes  in  Figure  5.3  are  root  mean  square  values  and  that  the 
calculated  data  points  in  Figure  5.3  refer  to  a  single  value  of 
0/ aQ.  [For  different  values  of  a/aQ,  one  might  naturally  expect  to 
obtain  a  family  of  curves  similar  to  Figure  5.3]. 


It  is  observed  from  Figure  5.3  that  the  total  entrainment  rate  at 
any  oscillation  frequency  always  lags  the  foil  motion.  The  maximum 

phase  angle  4^,  is  seen  to  be  about  20  degrees.  For  large  values 

of  k',  the  amplitude  of  the  nonsteady  entrainment  rate  tends  to 

decrease  rapidly,  as  shown  by  the  dotted  line  in  Figure  5.3,  with  the 

limiting  value  of  zero  for  very  high  frequencies.  At  least,  the 

present  experimental  data  for  reduced  frequencies  k'  up  to  k*  ■  1.0 
would  seem  to  indicate  such  a  trend.  The  rapid  decrease  in  cavity 
length  amplitude  at  reduced  frequencies  k'  higher  than  0.23  in 
conjunction  with  the  time  resolution  available  from  the  video-tape 
records  used  for  measuring  cavity  lengths  made  it  impossible  to  obtain 
cavity  length  amplitudes  for  frequencies  higher  than  k*  •  0.23.  It 
is  of  interest  here  to  observe  that  the  nonsteady  entrainment  rate 
lags  the  foil  oscillation  by  nearly  90°  at  all  frequencies.  It 
appears  that  as  the  cavity  is  increasing  in  volume,  the  cavity  is  able 
to  sustain  higher  cavity  pressures  and  thus  lower  entrainment  results 
behind  the  cavity.  The  uncertainty  involved  in  obtaining  data  points 
in  Figure  5.3  is  difficult  to  predict,  given  the  many  variables  and 
assumptions  involved  in  the  calculations.  However,  the  trend 
predicted  by  the  semi -empirical  model  is  interesting  and  should  serve 
as  a  guide  for  future  research  in  this  area. 

5.5  Recapitulation 

In  this  chapter,  the  two-dimensional  gaseous  cavity  closure  in 
steady  as  well  as  nonsteady  cavity  flow  is  analyzed.  It  is  observed 
that  the  presence  of  gases  within  the  cavity  alters  the  gas 
entrainment  rate  behind  the  cavity.  The  change  in  steady  gas 


entrainment,  however,  is  seen  to  be  insigifnicant  as  the  ratio  a/a 
decreases.  By  imposing  the  fact  that  the  surface  of  a  nonsteady 
cavity  is  essentially  a  material  surface  at  every  instant  and  assuming 
that  the  surface  waves  from  the  upper  surface  and  lower  surface  of  the 
cavity  terminate  at  one  point  behind  the  cavity  at  every  instant,  the 
phase  relation  between  the  foil  motion  and  cavity  motion  is  obtained. 
This  phase  relation  is  seen  to  be  consistent  with  the  experimental 
observations.  The  nonsteady  gaseous  entrainment  behind  the  cavity  is 
estimated  using  the  nonsteady  continuity  equation,  with  the 
assumption  that  the  gas  within  the  cavity  undergoes  isothermal 
compression  and  expansion.  The  nonsteady  gaseous  entrainment  is  seen 
to  depend  primarily  on  the  cavity  volume  fluctuations  and  cavity  gas 
pressure  fluctuations  because  it  is  found  that  for  the  range  of 
reduced  frequencies  considered  here  the  diffusion  into  the  cavity  is 
quasi-steady.  It  is  also  seen  that  the  nonsteady  entrainment  rate 
always  lags  the  foil  oscillation. 
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CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  the  aim  of  this  thesis  to  address  and  resolve  some 
of  the  aspects  involving  gaseous  diffusion  in  steady  and  unsteady 
cavity  flows  so  that  the  phenomenon  of  gaseous  entrainment  behind 
a  cavity  may  be  better  understood  and  an  improved  representation  for 
the  influence  of  the  nonsteady  cavity  in  a  water  tunnel  test  section 
can  be  obtained. 

We  have  demonstrated  very  clearly  that  existing  mathematical 
models  for  gaseous  diffusion  across  axisymmetric  cavity  surfaces  are 
inadeauate.  While  the  exact  solution  of  the  three  coupled  integral 
equations  formulated  in  Chapter  2  might  result  in  a  more  accurate 
representation  of  the  gas  diffusion  phenomenon  across  axisymmetric 
cavity  surfaces  than  has  been  possible  until  now,  it  seems  likely 
that  a  more  realistic  model  of  the  important  physics  would  be  more 
fruitful.  Even  so,  the  highly  simplified  approach  used  with  the 
best  experimental  results  available  from  the  literature  appears  to 
give  satisfactory  estimates  of  diffusive  mass  flows  into  the  cavity 
from  the  free  stream. 

As  noted  above,  a  plausible  alternative  for  the  solution  of  the 
axisymmetric  diffusion  problem  is  by  a  two-layer  diffusion  model. 
Such  an  asymptotic  solution  of  the  problem  will  represent  a  more 
realistic  physical  model  of  the  phenomena  which  could  account  for 
molecular  gaseous  diffusion  from  the  free  stream  into  the  turbulent 
concentration  layer  and  turbulent  diffusion  from  the  turbulent 
momentum  layer  into  the  cavity. 
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From  the  experimental  results  obtained  at  upstream  and 
downstream  locations  with  respect  to  the  unsteady  cavities,  we  have 
been  able  to  Identify  and  quantify  the  unsteady  cavity  as  a  dynamic 
element.  The  cavity  behind  an  oscillating  foil  has  been 
characterized  as  a  dipole  source  at  low  frequencies  of  oscillation. 
Unsteady  cavity  length  measurements  have  shown  that  for  sinusoidal 
oscillations  of  the  body,  higher  harmonics  of  the  cavity 
oscillations  are  important  to  be  considered,  especially  the  first 
two  harmonics  are  seen  to  be  significant.  Steady  cavity  pressures 
are  generally  seen  to  be  much  higher  than  the  vapor  pressure. 

Furt  *;r  experiments  are  definitely  warranted,  perhaps  in  a 
larger  water  tunnel,  to  extend  the  range  and  scope  of  the  present 
experimental  investigation.  The  measurement  of  phase  difference 
between  body  motion  and  cavity  motion  for  several  reduced 
frequencies  would  be  most  helpful  in  modeling  unsteady  cavity 
closure.  Even  so,  the  unsteady  gas  diffusion  analysis  presented  in 
Chapter  4  in  conjunction  with  the  experimental  results  has  enabled  a 
more  realistic  cavity  closure  model  to  be  developed  than  has  been 
possible  heretofore. 

The  steady  cavity  dossure  for  gaseous  cavities  has  shown  that 
the  presence  of  noncondensable  gas  in  the  cavity  does  affect  the 
steady  entrainment  rate  of  gases  at  the  cavity  terminus.  The 
formulation  of  nonsteady  gaseous  cavity  closure  involves  the 
determination  of  phase  angles  between  the  foil  motion  and  i) 
nonsteady  gas  diffusion  through  the  cavity  surface,  ii)  the  cavity 
volume  fluctuations,  and  iii)  cavity  pressure  fluctuations, 
respectively.  By  a  simple  cavity  terminus  closure  model,  we  have 


been  able  to  estimate  the  phase  difference  that  could  exist  between 
the  foil  motion  and  the  cavity  length  motion.  The  many  variables 
Involved  in  estimating  the  nonsteady  entrainment  rate  clearly 
suggests  the  need  for  more  refined  experimentation  and  theory  to 
understand  more  fully  nonsteady  gaseous  cavities  and  the  related 
turbulent  entrainment  rates. 
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A. 1  Asymptotic  Evaluation  of  Integral  in  Equation  2.10 
The  integral  in  Eq.  (2.10)  is 


2t  -x  V a"  -  2rcos8 


dd  , 


(A- 1.1) 


a  -  2rcos9 


where 


2  (  _  >  2  2  . 
m*(x-£J+r+l 


(A. 1.2) 


(We  recall  here  that  the  distances  r,  x  and  %  are  non-dimensionalized 
with  respect  to  the  radius  of  the  cavity  and  X  -  Um  a/ 2k  is  a  large 
number  of  the  order  1000.) 

We  use  Laplace's  method  (see,  for  example,  Nayfeh,  A.  H. , 
Introduction  to  Perturbation  Techniques,  John  Wiley,  1981,  p.  88) 
to  evaluate  the  integral  in  Eq.  (A. 1.1). 

Letting 


(A. 1.3) 


m  -  2rcos0 


h(8)  “  -  V®  “  2rcos0  . 

The  integral  in  Eq.  (A. 1.1)  may  be  written  as 
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I  -  /  f(0)eXh(9)d0  . 
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We  will  now  expand  h( 0)  and  f ( 9 )  in  Taylor's  series, 
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so  that  ^—7} 
d9~ 


9-0 


Hence  the  function  f(9)  also  has  an  absolute  aaxiauQ  value  at 
9-0.  f(9)  may  now  be  expanded  in  a  Taylor's  series. 


f  (  9) 


r9 


2(a2  -  2r)3/2 


(A. 1.7) 


Substituting  the  series  expansion  for  f ( 9 )  and  h(9)  in  Eq.  (A.1.5), 
we  have 
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(A. 1.8) 


In  Eq.  (A.  1.8),  the  Halt  -  instead  of  lit  appears  because  the 
absolute  aaxiaua  value  of  the  function  h(9)  occurs  at  9  -  0  and 
hence  the  aaxiaua  contribution  to  the  integral  results  froa  the 
region  in  the  neighborhood  of  9-0.  Equation  (A. 1.8)  Bay  be 


rearranged  as 
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A. 2  Simplification  of  Integral  Equations 

Simplification  of  Eqs.  (2.21),  (2.22)  and  (2.23):  Consider 
Eq.  (2.21).  The  third  integral  on  the  right-hand  side  of  Eq.  (2.21) 


-21(5  -  x  ) 

1  f,U)e 

/  - - - - d?  . 

*2  “  *2 

Recognizing  the  fact  that  X  is  a  large  number,  the  strong 
negative  exponential  causes  the  function  f^  to  be  a  slowly  varying 
function  and  thus  f^  may  be  assumed  to  contribute  to  the  integral 
only  near  5  ■  x^.  Thus, 


1  f,(5)e 


-2X(?  -  x  } 


^5  - 


£,  “ 2X(5  — 
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The  upper  limit  l  may  be  replaced  by  •  and  the  integral 
evaluated  to  get  ^(*2  The  very  same  arguments  lead  to  the 

simplification  of  second  integral  of  Eq.  (2.22)  and  fourth  integral 


of  Eq.  (2.23). 


Similarly,  the  fourth  integral  on  the  right-hand  side  of 


Eq.  (2.21)  is 
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Letting  5  -  -  a,  we  have 


Integral  ■  /  - — - da  . 

»_»  /a 


Once  again  because  of  the  strong  negative  exponential,  we  nay  write. 
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We  now  note  that  is  defined  only  in  the  interval  £  <  5  <  •. 


Thus,  the  fourth  integral  on  the  right-hand  side  of  Eq.  (2.21)  is 


zero.  The  fourth  integral  on  the  right-hand  side  of  Eq.  (2.22)  is 
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By  leeting  5  -  x  -  a,  we  have 
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A. 3  Solution  of  Eo .  (2.22a) 
Equation  (2.22a)  is 


/  -  d*  +  f  i^xi^  T 


(A.3.1) 


We  rewrite  Eq.  (A.3.1)  as 
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where 
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We  shall  use  Picard's  method  with 
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(A.3.4) 


where  p  and  c^  are  unknown  constants.  Then  substituting  Eq.  (A.3.4) 
in  Eq.  (A.3.2), 
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If  we  now  let  t  *  -  5,  we  have 
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We  note  Chat  j(x^)  i*  not  *ery  different  from  f^  0(x^). 
we  choose  p  such  that  g(x^)  *  f^  ^us  w*  obtain 

."•ci  *  •*  v?'  ■  °  • 

Talcing  the  negative  square  root*  so  p  is  real,  we  get 

i"‘ V^’° 


or 
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x^  <  0  . 


Hence 
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*This  is  required  for  the  function  fj(xi)  to  satisfy  the  integral 
equation. 


A. 4  Solution  of  Eg.  (2.21a) 


Equation  (2.21a)  is 
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Substituting  the  value  of  from  Eq.  (2.24),  we  have 
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Consistent  with  our  approximations  in  the  first  integral  appearing 
on  the  right  hand  side  of  Eq.  (A.4.2),  we  shall  assume  that  the 


strong  negative  exponential  will  cause 


1 


to  vary  slowly  so 


/x2  -  C 

that  it  can  be  assumed  to  contribute  to  the  result  only  near  £  -  0. 
Thus, 
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Taking  the  Laplace  transform  and  simplifying,  we  have 
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where  F2  is  the  Laplace  transform  of  f^.  Taking  the  inverse 


Laplace  transform,  we  find  for  f^, 
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APPENDIX  B 


EXPERIMENTAL  RESULTS  AND  ANALYSIS 


B.2  Measurement  of  Free  Air  Content  in  the  Test  Section 


The  change  in  acoustic  velocity  as  a  function  of  free  gas 
concentration  is  given  as  (see  for  example,  John  M.  Killen  and 
John  F.  Ripken,  "A  Water  Tunnel  Air  Content  Meter,"  St.  Anthony 
Falls  Hydraulic  Laboratory  Report  No.  70) 


where 


(3.2.1) 


C %  *  velocity  of  sound  in  the  mixture 

Cy  *  velocity  of  sound  in  water 

x  ■  concentration  of  free  gas  by  volume 

Ew  *  bulk  modulus  of  pure  water 

Pa  «  absolute  pressure  in  the  test  section. 


In  order  that  Eq.  (B.2.1)  applies,  it  must  be  established  that  the 
gas  bubble  sizes  present  in  the  test  section  are  smaller  than  the 
resonant  size  corresponding  to  the  frequency  of  the  measuring 
acoustic  signal.  From  Eq.  (B.2.1)  we  can  obtain  an  expression  for 
the  concentration  of  gas, 


x 


c2 

,(4  -  o 


(B.2. 2) 


B.3  Analysis  of  Experimental  Error 


( 1)  Cavitation  Number 

The  cavitation  number  a  is  defined  as 
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n  * 
l/2o  U“ 

oo  cm 


where 


(B.3. 1) 


p0  *  test  section  pressure 
pc  ■  steady  cavitation  pressure 
*  density  of  tunnel  water 
U«  ■  free  stream  velocity  in  the  test  section. 


From  Eq.  (B.3. 1) , 
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(B.3. 2) 


where  the  A  before  a  quantity  signifies  a  small  error  involved  in 
measuring/estimating  that  quantity. 

The  maximum  error  in  calculating  (p0  -  pc)  is  0.12  psi.  The 
minimum  value  of  (pQ  -  pc)  during  the  experiments  is  seen  to  be 
2.7  psi.  The  maximum  error  in  calculating  the  free  stream  velocity 
is  0.05  ft/sec.  The  minimum  value  of  U*  during  the  experiments 
is  25  ft/sec.  With  these  extreme  values,  the  maximum  absolute 
error  in  estimating  the  cavitation  number  may  be  calculated  as 


Hence  the  error  la  estimating  the  cavitation  nunber  is  at  the  most 


(2)  Concentration  of  Sea 

From  Eq .  (B.2.2)  the  concentration  of  gas  is  given  by 
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t 


where 


Pa  ■  test  section  pressure 

*  bulk  modulus  of  pure  water 
Cv  “  speed  of  sound  In  pure  water 
Cg  *  speed  of  sound  In  the  mixture. 
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We  want  to  bound  the  error  In  the  estimation  of  concentration  of 
gas  when  there  is  a  known  uncertainity  in  the  measurement  of  speed 
of  sound  in  the  mixture. 

From  Eq.  (B.2.2), 
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noting  that  Ew  »  Pa 


The  maximum  error  in  measuring  Che  speed  of  sound  in  Che 


mixture  is  30  £c/sec.  The  minimum  soeed  of  sound  in  Che  mixture 


during  Che  experiaencs  is  seen  Co  be  abouc  4000  fc/sec.  The  error 


in  reading  the  tunnel  pressure  is  0.1  psi  and  Che  minimum  tunnel 


pressure  is  3  psia.  Hence  from  Ea.  (3.3.3), 
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Therefore  the  bound  on  the  error  in  estimating  Che  concentration 


of  gas  is  8%. 
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B.4  Wall  Effect  Correction  for  Cavitation  dumber 

The  technique  described  by  Hsu  (24)  is  used  in  obtaining  the 
true  cavitation  number.  Following  Hsu,  the  linearized  results  for  a 
two-dimensional  flat  plats  hydrofoil  with  trailing  cavity  in  a 
closed  tunnel  are: 


(B.4.1) 


where 
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to 
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(B.4. 2) 
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In  the  above  equations. 


a  *  angle  of  attack 
c  ■  chord  length  of  foil 
l  »  cavity  length 


L  -  l  +  c 


H  ■  tunnel  diameter. 
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(B.4.3) 


Having  found  <jj  from  Eq.  (B.4.1),  the  true  cavitation  number 
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is  found  from 
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Example :  The  cavi cation  number  calculated  from  the  measured  cavity 
pressure  is  0.75.  The  cavity  length  l  is  7.32  cas.  The  chord 
length  of  the  foil  is  7.14  cas.  The  angle  of  attack  is  10  degrees. 
From  So.  (B.4.2), 


h  '  2(l§o  “HthI  -  *1 


Xf  -  0.353  . 


From  Eq.  (B.4.3)  for  *  2.981, 
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w  Z 


\  -  2.34  . 
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From  Eq.  (B.4. 1) , 


•i  -  2-34  *  1-277 


From  Eq.  (B.4.4) , 
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C. 1  Fundamental  Solution  of  So .  (4.13b) 

Fundamental  solution  to  Equation  (4.13b)  for  an  oscillating  unit  source 
at  the  origin: 


Equation  (4.13b)  is 
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(C.1.1) 


We  assume  that  the  cavity  length  fluctuations  are  harmonic.  Then,  the 
entire  concentration  field  Cy  is  harmonic  so  that 


cu(x,y,t)  -  c  (x,y)e- 


(C.1.2) 


where  m  is  the  circular  frequency  of  harmonic  oscillations  and  j  =*  /-I  .  It 
is  to  be  noted  that  only  the  real  part  of  the  analysis  that  follows  is 
physically  pertinent.  Substituting  Equation  (C.1.2)  in  (C.1.1),  we  obtain 
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We  normalize  the  x  and  y  coordinates  by  the  flat  place  chord  b  and  simplify 


Eq.  (C.1.3)  as 
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7  -  2-  ,  (C.1.6) 

b 


k  »  ~  m  reduced  frequency,  (C.1.7) 
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R  ■  — —  »  Diffusive  Sevnolds '  number  (C.1.3) 
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To  solve  Equation  (C.1.4)  we  replace  the  Laplacian  on  the  right-hand  side  by 
its  counterpart  in  plane  polar  coordinates  when  variations  only  in  the  radial 
direction  are  permitted.  Then, 
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(C.1.10) 


It  oust  be  noted  that  the  coordinates  x  and  r  appearing  in  Eq.  (C.1.9)  are 
dimensionless. 

By  letting 


(C.1.11) 


Equation  (C.1.9)  can  be  written  as 


(C.1.12) 


We  now  use  separation  of  variables  to  solve  Equation  (C.1.12). 
substituting 

C(s,r)  -  S(s )X(r )  , 

Equation  (C.1.12)  aav  be  written  as 


Bv 


(C.1.13) 


t  V  f  V 

...  „  ^  S  X  ^  1  X  2 
2jk  R  +  —  -  —  +  -  < 


X 


1  ’ 


(C.1.14) 


where  is  an  unknown  constant.  Therefore, 


*  *  x  2 

X  +  —  -  tef  X  -  0  , 


(C.1.15) 


t  t 

S  2 


f—  -  <1  -  2jkR 


(C.1.16) 


Equation  (C.1.15)  may  be  rewritten  as 


2  d  X  dX  2  L  . 
r  _  +  r  -  -  .  -  0  , 

dr 


(C.1.17) 


and  the  solution  for  ic^  >  0  is 


X(r)  -  A,  I  (k  r)  +  A  K  (ic.r)  ,  (C.1.18) 

10  1  £01 

where  A^  and  A2  are  constants  and  Kg  and  I0  are  the  modified  Bessel's  function 

of  zeroth  order  of  first  and  second  kind  respectively. 

2 

Note:  Negative  values  of  ic^  yield  functions 


e-r_ 
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Jo  and  Y0  (Bessel's  function  of  zeroth  order  of  first  and  second  kind 
respectively)  which  do  not  satisfy  the  boundary  conditions. 

2 

Uote  2:  ■  0  yields  a  logarithmic  function  which  again  does  not  satisfy 

the  boundary  conditions. 

For  large  distances  from  the  origin,  the  function  c(s,r)  should  vanish. 
Hence,  we  sec  Aj  *  0  in  Equation  (C.1.18).  Therefore, 


X(r)  -  A,K.(<tr). 

Z  o  1 

The  solution  of  Equation  (C.1.16)  is  easily  written  as 


(ie‘  -  2jkR)s 


S(s)  *  Be 


where  B  is  an  unknown  constant.  Hence,  from  Equation  (C.1.13), 


(C.1.19) 


(C. 1.20) 


I 


(<  -  2jkR)s 

C(s ,r)  -  Bl  e  Ko(ielr)  * 


(C.1.21) 


where  Bj_  *  BAo  is  a  new  unknown  constant.  The  arbitrary  constant  and  the 
separation  constant  are  obtained  by  requiring  that  Eq.  (C.1.21)  go  to  the 
proper  steady  state  limit  (8)  when  the  reduced  frequency  k  goes  to  zero.  This 
results  in 


B 


and 


1  2rpv 


■  R 


Hence,  for  a  unit  oscillating  source, 


C<X’C)  "  22^  e"jkX  ^  Ko  (Rr) 


(C. 1.22) 


(C.1.23) 


(C.1.24) 


C.2  Inverse  Laplace  Transfora 

To  obtain  the  inverse  Laolace  Transfora  of  the  function 


Denote 


F(s,t)  - 


_ 0  /s  +  jk _ 

s  [  /2R  -  jk  +■  /s  +  jk] 


a  »  jk, 


b  »  /2R  -  jk  . 


F(s,t) 


0  /s  +  a 

s  [/rrr  +  b] 


F(s,t)  ,  Q[JSL+-*+  b - - ]  , 

s  [/s  +  a  +  b]  s[/s  +  a  +  b] 


P(i.t')  ■  Q(l-  _J_ - -]  . 

s [/s  +  a  +•  b] 


8  [/s  +  a  +  b] 


b  f  Vs  +  a  —  b 


s  [s  +  a  -  b  ] 


a  ■  a  -  b  . 


(C.2. 1) 


(C.2. 2) 


(C.2. 3) 


(C.2. 4) 


(C.2. 5) 


(C.2. 6) 


(C.2. 7) 


(C.2. 8) 


?#*£ 

Jii? 


i 

m 


P 

1; 


|*w,* 

wi 

5?: 


s(s  +  o)  /s  +  a 


s  /s  +  a  s(s  ■+■  a)  /s  +  a 


3  /s  +  a  o  s  /s  +  a  (s  +  a)  /s  +•  a 


Substituting  Equation  (C.2.11)  in  (C.2.9),  one  has 


s  /s  +■  a  as  /s  +  a  a  (s  +  a)  /s  +  a  s(s  +  a) 


(C.2.12) 


Substituting  Equation  (C.2.12)  in  Equation  (C.2.6),  one  gets 


F(s,t)  -  Q[i  - 


s  /s  +  a  as  /s  +  a  a(s  +  a)  /s  +  a  s(s  +  a) 


(C.2.13) 


Talcing  the  inverse  Laplace  Transform  on  both  the  sides,  one  finds 


f(x,t)  ■  Q[—  -  —  e  ^erfcb  /x  -  erf  /ax]  . 

a  a  a 


(C.2.14) 


«y 

$;;l 

m 


m\ 

ml 


CVS' 

BW 

c»*  y 
EW»» 

|.SV*  I 


rly,/*  I 
[■$$ 


Mfl 

raVM.V.M 


C.3  Integral  Evaluation 


,V| 

m 


I  >'J 

I  >3 


hi 

I'll 


Equation  (4.25)  is 


i(t) 

/  f(x,t)dt  , 
0 


(C.3. 1) 


where 


f(x,t)  *  —  (a  -  b~e  ^erfcb  /x  -  b  /a  erf  /ax) 


(C.3. 2) 


Substituting  Equation  (C.3. 2)  in  (C.3. I), 


i(t)  £(t)  _  _i(t)  __ 

-  —  [  /  adx  -  b ^  /  e  ^erfcb  /x  dx  -  b  /a  /  erf  /ax  dx] 
' z  a  0  0  0 


(C.3. 3) 


t(t) 

/  e  ^erfcb  /x  dx  , 
0 


(C.3.4) 


£(t) 

»  J  erf  /ax  dx  . 
0 


(C.3. 5) 


The  integral  in  Equation  (C.3.4)  may  be  evaluated  by  parts. 


-ox  i(t)  ,  -b  x 

I.  ■  erfcb  /x  - - ]  +  —  /  e-aX(-l)  - — - dx 

1  -  a  1  n  a  ‘  r 

0  0  /it  /x 


m 


m 


m 


0 

v,;.* 


m 

Ml 


By  setting  ax  »  n^  in  Equation  (C.3.6),  we  can  show  that 


_1_  e  a^arfcb  /£ 
i  "  a  '  a 


erf  /ai  . 


a  /a 


The  integral  in  Equation  (C.3.5)  can  be  obtained  by  parts. 


£(t)  *(0  jT  ^~ax 

1  -  xerf  /ax]  -  /  x - dx 

2  0  n  /»  /x 


(C.3.7) 


1  *6  •  | 


*2 


ierf  /a  l  - 


-  lit) 

~h{ 


-ax, 

/x  e  dx 


(C.3.8) 


The  integral  in  Equation  (C.3.8)  may  be  evaluated  again  by  parts  and  the 
result  is 


I.  *  £erf  /a l  +  —  e  a*  -  -I-  erf  /a£ 
2  na  2a 


(C.3.9) 


Hence  from  Equation  (C.3.3), 


g.  .  a  (a4  .  i_  +  i_  e-»»er£cb  fT  *  -S—  erf  /  ai 


dt 


a  a 


/r 


2  /  a 


erf  /  a£  -  b  —  e-3*  -  b  /  a  £erf  /  a£  ] 

IT 


(C.3.10) 


APPENDIX  D 


SENSITIVITY  ANALYSIS  FOR  MASS  DIFFUSION  RATE 


'l*  4* 


I  '  |  •  «L*  *  ' 

i  *  ^  *  t. 


to; 

SsS- 


I 

i 


®,i 

$ 

IW. 


tea 


I*'*.* 


K® 


B® 

Sw 

yw 

m\ 


Kj*I 

M*. 

ft!* 

P'M 


The  quasi  steady  mass  flow  rate  across  a  two-dimensional  cavity 


is,  from  Equation  (4.23) 


mqs  ■  iAI  *  c0  c 


2R2,  - 

e  errc 


2R 


VirT 

Y  *R  ~  2x 


(4.28 


For  large  values  of  the  argument  2R£, 


i2i  erf c  2RZ  »  —  ( 
/it  ' 


.  3/2 


,  +  ...  )  .  (D.l) 


2R£  2(2Ri) 


Substituting  Eq.  (D.l)  in  Eq.  (4.28)  above. 


4AT  co 

/  TT 


(V2R  l  -  -2  + 


•  — 1 —  3/,  +  ...)  .  (d.2: 

2V2RZ  4(2R2) 


Retaining  only  the  first  term  within  the  brackets  in  Eq.  (D.2),  we 
take  the  logarithm  on  both  sides  and  then  differentiate  to  obtain 


4c 

o 


c 

o 


(D.3) 


and  this  is  the  basis  of  the  order  of  magnitude  discussion  surrounding 


Eq.  (4.31)  of  Chapter  4 


APPENDIX  E 


UNSTEADY  ENTRAINMENT  DATA  AND  CALCULATIONS 
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E.2  Quasi-Steady  Concentration  Fluctuation  Amplitude 
Equation  (5.2)  is 


-  0.0976  —  +  0.74 
a 


(5.2) 


Taking  logarithms  on  both  sides  and  then  differentiating  Eq.  (5.2), 
one  obtains 


Ap„  ■  0.0976(a/a  )(Aa  /a  )  p  .  (E.2.1) 

*  G  0  0  0  V 

o 

From  Eq.  (2.9),  the  concentration  difference  c  expressed  in  ppm  by 
modes  is 


so  that 


9 


apg  /b 

o 


c 

o 


Substituting  from  Equation  (E.2.1), 


Ac  -0.0976  (a/a  )(Aa  /a  )(P  / B) 

o  m  _ o _ oo  v _ 

c  a,  -  (p  /S)(0.74  -  0.0976  a/a  ) 

O  1  V  O 


(2.9) 


(E.2. 2) 


E.3  Evaluation  of  3m/ 3t  in  Eg.  (5.16b) 


Let  V  ■  total  volume  of  cavity  in  cubic  feet, 

m  »  total  mass  of  gas  in  the  cavity,  in  slugs. 

Then,  at  any  instant,  the  equation  of  state  for  the  gas  within  the 


cavity  is 


pQV  -  mRT 


(E.3. 1) 


Where  R  is  the  specific  gas  constant  and  R  *  1718  ft. lb. /slug  °R 
for  air. 

Assuming  that  the  temperature  of  gas  within  the  cavity  remains 
constant,  Eq.  (E.3.1)  may  be  partially  differentiated  with  respect  to 
time  to  obtain 


3m  1  ^  3 V. 

3t  "  RT  V3t  PG  3t^ 


From  Eq.  (5.13a),  the  cavity  volume  V  is 


v,  if i  (I -i) 


where  a  ■  a  (1  +  ActeJ  ) 
o 


(E.3. 2) 


(5.13a) 


i  -  I  (1  +  ceJ(“C++))  . 


Differentiating  Eq.  (5.13a)  partially  with  respect  to  t  and 
neglecting  terms  of  order  eAa,  one  obtains  after  suppressing  the 


e^ut  factor 


<* *bk'qJ  -  u  i 
- = -  [i  2  +  (t  -  4-)  Aa]  . 


(E.3. 3) 


»Y *%.»'»  4 


Now,  from  Eq.  (A. 17),  the  cavity  gas  pressure  pG(t)  is 


j  (wt+<p  ) 

Pg(t)  -  PG  (1  +  5e  a  )  .  (4.17) 

o 

Differentiating  Eq.  (4.17)  partially  with  respect  to  t  and 
suppressing  e^^  factor. 


an  P r  u 

3pG  Go  C  j*a 


at 


(E.3.4. ) 


Substituting  Eqs.  (E.3.3)  and  (E.3.4)  in  Eq.  (E.3.2),  one  obtains 


imb  pQ  k’qcj  J+  _ 

- - -z— -  [U-^KSe  a  +  Aa)  +  l  . 


(S.3.5) 
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E.4  Unsteady  Entrainment  Calculated  Data 

The  sample  calculations  below  show  how  the  tabulated  results  were 
obtained  for  one  case  at  a  10  Hz  frequency  and  2°  peak  to  peak  foil 
amplidtude. 

From  Eq.  (5.3) , 


4pv  c 


d  j— 

O  1  IT 


—  [Vzr i  +  V:rI,  ]  . 


(5.3) 


R  -  yy  *  2500  ,  c  -  10  ppm,  o/aQ  ”  2.7, 

L  «  2.32,  -  1.32,  M.  -  1.46  x  10-5  slugs/sec  ft 

u  l  d 

o  00 


From  Eq.  (E.2.2) , 


l 

u  l 
o  o 


-  -0.0025 


0.165  x  0.707 


2  /  2.32  x  1.32 


0.033 


The  last  term  in  Eq.  (5.21)  is 


..0bP(.k'qc3  „  . 

- 2 -  [U,  -  ±  )(ia  +  Se  a)  +  t  re3*] 

2RT  A.  "o  2  “o 

d 

o 

From  Eq.  (5.5)  $  -  k*  ■  0.03  radians  ■  1.72°.  From  Equation  (5.2) 

and  (E.2.1),  6  -  0.05. 


ft* 

tv  « 

‘vM 


ft: 

KW 


4 

m 

I'V 

y,y 

■r 


% 

1 

i(r| 

v.v 

5# 


144 


From  Eq.  (5.15) , 

Ao  . 

v  -  (ee^  +  Act - =—  )/(l - z - ) 

21  21 

u  u 

O  0 

v  -  0.239  +  0.006j  -  0.239  @  1.44° 

$a  *  180°  +  1.44  =*  181.44°  . 

From  Fig.  (3.14),  for  a! aQ  ■  2.7,  *  72  lb/sq.ft. 

K'  o 

Substituting  the  values,  the  last  term  in  Eq.  (5.21)  is 
-0.002  +  0.065j,  as  tabulated. 
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